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Abstract 



We continue our study of the decoupled wave equation in the exterior of a spherically symmetric, 
Schwarzschild, black hole. Because null geodesies on the photon sphere orbit the black hole, extra effort 
must be made to show that the high angular momentum components of a solution decay sufficiently 
fast, particularly for low regularity initial data. Previous results are rapid decay for regular (H :i ) initial 
data [6] and slower decay for rough (H 1+€ ) initial data [5]. Here, we combine those methods to show 
boundedness of the conformal charge. From this, we conclude that there are bounds for global in time, 
space-time norms, in particular 



for H +e initial data with additional decay towards infinite and the bifurcation sphere. Here <j> refers to 
a solution of the wave equation. I denotes the exterior region of the Schwarzschild solution, which can 
be expressed in coordinates as r > 2M, t £ t,ii) £ S 2 , and d 4 vol is the natural 4-dimensional volume 
induced by the Schwarzschild pseudo-metric. 

We also demonstrate that the photon sphere has the same influence on the wave equation as a 
closed geodesic has on the wave equation on a Riemannian manifold. We demonstrate this similarity 
by extending our techniques to the wave equation on a class of Riemannian manifolds. Under further 
assumptions, the space-time estimates are sufficient to prove global bounds for small data, nonlinear 
wave equations on a class of Riemannian manifolds with closed geodesies. We must use global, space- 
time integral estimates since L°° estimates cannot hold at this level of regularity. 

1 Introduction 

The Schwarzschild solution is a manifold which satisfies the Einstein equation. Recall that, in general rela- 
tivity, space-time is described by a four dimensional, Lorentz manifold which satisfies the Einstein equation. 
In vacuum, this equation reduces to the statement that the Ricci curvature is zero. 

In 1918, Schwarzschild found the first non-trivial solution, which is the unique spherically-symmetric, 
vacuum solution, other than Minkowski space, flat R 3+1 with the metric — dt 2 + dx 2 . The Lorentz pseudo- 
metric for the Schwarzschild solution is 



This is the pseudo-metric outside any spherically symmetric body with mass M, and, for large r, the space- 
like geodesic motion on this manifold, as described by the coordinates (t, r, 9, </>), approaches that of a body 
moving under the influence of a central mass, according to Newton's laws, as described by polar coordinates in 
R 3+1 . Similarly, null geodesies describe the motion of light rays, which physicists treat as massless particles, 
called photons. 

The pscudo- metric (|1.1[) is singular at r — and appears to be singular at r = 2M, and this was 
believed to be the case for many years; however, this is an artefact of the coordinate system, and in the 
completion of the manifold, the surfaces corresponding to r = 2M are null geodesic surfaces. The structure 




ds 2 = -(1 - 2M/r)dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 {d9 2 + sin 2 (6>)^ 2 ) 



(1.1) 
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Figure 1: The causal structure of the Schwarzschild solution. 

of the completion of the Schwarzschild solution is a little surprising, but, by now, well understood [l4j, and 
illustrated, up to a conformal transform and ignoring the angular coordinates, in figure [I] 

Since as r — * oo, the Schwarzschild solution approaches Minkowski space, the causal boundary near 
r — > oo is like that in Minkowski space, with a future and past time like infinity, a spatial infinity, io; 
and future and past null infinities, 3 ± , which corresponds to a null hyper-surface when the Schwarzschild 
solution is embedded via a conformal transform (ds 2 i— ► VL 2 ds 2 ) in a compact Lorentz space. The intersection 
of the chronological past of the future null infinity with the chronological future of the past null infinity, 
P\I + (3~) is commonly referred to as the exterior region of the Schwarzschild black hole and depicted 
as region I in figure [TJ The boundaries of the past of the future null infinity and of the future of 

the past null infinity both correspond to surfaces where r — 2M. The exterior region is described 

by the Lorcntzian pseudo-metric (jl.ip on the set 

I = R x Tl = {(t,r,9,4>)\t £R,r> 2M, {9,<p) £ S 2 } 

The singularity at r = corresponds to both the space like surfaces S + and S~ in figure [1] The regions 

II and IV around them are both described by the Lorentzian pseudo-metric fll.lj) with r < 2M. Region 

III is completely isomorphic to region I. The surfaces r = 2M are null hyper-surfaces around the region 
r < 2M which are referred to as the black hole, since light rays can enter but not escape from them. More 
precisely, since events inside the black hole (in region II) can not be witnessed by observers remaining outside 
(headed to i + or the surfaces r = 2M are referred to as event horizons. The coordinate t diverges on 
the event horizon (diverging to negative infinity in the past, and plus infinity in the future), so that it fails 
to distinguish points on the event horizon. The points where the four regions meet, corresponding to taking 
the limit, in region I, as r — > 2M for any fixed value of t, is called the bifurcation sphere. 

At r = 3M, there are null geodesies which orbit the black hole, remaining at r = 3M. Since photons 
can orbit here, this surface is known as the photon sphere. Null geodesies initially tangent to a surface of 
constant r will escape to infinity if r > 3M and fall into the event horizon, and then into the singularity, 
if r < 3AI . From this perspective, r = 3M is a natural boundary. However, this surface does not play a 
special role in terms of causal geometry. 

A major open problem in relativity is the stability of black holes under small perturbations. This is 
important for understanding the behavior, interaction, and dynamics of black holes. Exact solutions exist 
for black holes with electromagnetic charge and angular momentum, as well as mass. Charged solution are 
given by the Reissner- Nordstr0m solutions. Black holes with angular momentum rotate and are described 



2 



by Kerr solutions, or Kerr-Newman solutions if they are also charged. If the charge and angular momentum 
are not excessive, an event horizon and causal infinities still exist, and the geometry of region I is not greatly 
altered; although, the nature of the singularities changes dramatically and there are infinitely many additional 
exterior regions, like region III in the Schwarzschild solution. By the "no hair theorem" , stationary black 
holes arc entirely characterized by their mass, charge, and angular momentum, and it is generally believed 
that black holes approach one of the known exact solutions. From Birkhoff's theorem, it is known that 
there can be no spherically-symmetric, vacuum (or vacuum-Maxwell) perturbations [14j . Since the Einstein 
equation is nonlinear, it is expected to couple higher spherical harmonic perturbations so that higher spherical 
harmonic perturbations to the Schwarzschild solution will introduce rotation, and that, although the Kerr 
solutions are expected to be stable as a class, the Schwarzschild solution alone will not be. However, since 
the Kerr solution has no global time- like Killing vector (even just in the exterior), its perturbations are 
extremely hard to analyze. Perturbation equations for Schwarzschild and Kerr-Newman have been derived 
by Regge and Wheeler [27], Zirilli [35], and Teukolsky [33], and all of these suggest black holes are stable. 
Whiting has shown there are also no unstable modes [34] . 

The only rigorous proof of nonlinear stability is that of Christodoulou and Klainerman, which proves 
the stability of Minkowski space [9] (see also [M] for a simplified proof). All the perturbation equations 
derived and the analysis of Minkowski stability have involved second order wave equations, so that we hope 
that information on the decay of the linear wave equation will provide information for the stability of Kerr 
solutions. 

The rigorous analysis of the wave equation on the Schwarzschild solution began with [12) which proved 
the L°° boundedness of solutions with initial data in H 2 . Scattering results have been proven for wave 
equations on the Schwarzschild p] and Kerr [TH] solutions. Decay for various types of fields, including the 
wave equation, around a Kerr black hole have been proven using spectral techniques [17j . A similar approach 
to ours has been used to prove decay for the wave equation on the Schwarzschild solution, including estimates 
along the event horizon, known as Price's lawpj]. Estimates along the event horizon are important not only 
for perturbation in region I, but also for understanding the nature of the singularity [10] 

We consider the decoupled wave equation on the exterior of the Schwarzschild black hole, 

□0=0, 0(0) =0 Ol 8t$(0)=fa, (1.2) 

where the d'Alembertian □ is given by the fixed background pseudo- metric Ijl.ip . By introducing the Regge- 
Wheeler tortoise coordinate, r„, and transforming to = rcf> (and similarly for the initial data), we transform 
the wave equation to 

-d 2 = -d 2 t + F0+^(-A s2 )0, 0(0) =0o, 9*^(0) =0i. (1.3) 

This has the form of a family of one dimensional wave equations on each spherical harmonic, and we have 
previously treated it as such [S][S]. Our method is to introduce the energy and conformal charge, 

E=[ |9 t 0| 2 + l/|0| 2 +14|V S 20| 2 dr,d 2 w, (1.4) 

J{t}xWl 

E c = [ (t 2 +r 2 )(|a t 0| 2 + V^|0| 2 + F L |V 5 20| 2 ) + 2tr*(dt4>){d r „4>)dr*d 2 u + E. 

We analyze this using commutator methods and a radial derivative operator. Since the energy comes from 
a Killing vector, it has long been known to be conserved. The conformal charge has been used to analyze 
wave equations in R" +1 , for example [18]. Since the conformal charge, on the Schwarzschild solution, does 
not come from a Killing vector, it is not conserved. 

Previously, we have shown that the growth of the conformal charge can be estimated by the accumulation 
of the angular component of the energy which does not disperse from near the photon sphere, 

E c {h) - E c (t ) < / 1 f X (|0| 2 + \V S 20\ 2 )dnd 2 Lodt. (1.5) 

Here x is a compactly supported function coming from the positive part of the trapping terms for the 
potentials, 

2V + r*V',2V L +r*Vl. 
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Our analysis begins with a Morawetz estimat^ 

/ / X(I0| 2 + \d r «<P\ 2 )dr*d 2 Ludt < CE. (1.6) 
Jo Jm 

This estimate is proven using the technique first introduced in [25] but focuses on the role of different terms. 
Since angular derivatives can be applied to the linear wave equation (|1.3[) . we can apply the Morawetz 
estimate (|1.6[) to (| 1 . 5|) to get that the conformal charge is bounded in terms of (j> and the operator square 
root of (1 — As"), denoted L, by 

E c (t) <CtE[L<t>] +E C (0). (1.7) 

Since E is an L 2 norm of derivatives of </>, and L acts like an additional derivative V 52 , we refer to (| 1 . 7[) as 
an H 2 estimate. 

We are interested in improving (|1.7p since the growth in t of the conformal charge does not give good 
decay estimates for <f> or <j) an d the loss of angular regularity can make it hard to close nonlinear estimates. 
In [6] , it was possible to improve (|1.7p to uniform boundedness in time for the conformal charge by proving a 
Morawetz estimate localized inside the light cone and using t~ 2 Ec instead of E as the bound in the Morawetz 
estimate (|1.6|) . To extend to non- linear problems, L°° estimates were used, which we required an additional 
angular derivative, leaving a H 3 estimate. In [5], we were able to reduce the loss of angular regularity in (|1.6p 
from one derivative to e derivatives, which is defined by the spectral theorem. We did this by dilating the 
radial derivative operator by fractional powers of the angular derivative operator and by relating powers of 
the angular derivative to the radial derivative, techniques we refer to as angular modulation and phase space 
analysis. Here, we combine these techniques. Since we are working only at an H 1+£ level of regularity, it is 
not possible to prove an L°° estimate, and, instead we prove global-in-time, space-time integral estimates, 
the most transparent of which is 

Proposition 1.1. If <fi is a solution of (jl.2p and e > 0, then 

f \4>\ 4 d 4 vol < C{E C [0o, 0x] + E[^ , 0x] + E[L S <I> , L^]) 

J Region I 

provided the norms on the right are finite, the initial data vanishes sufficiently fast at the bifurcation sphere 
and spatial infinity, and where <p — r<p and d 4 vol is the four dimensional space-time volume associated to the 
pseudo-metric in p. ID . 

The natural measure d 4 vol can be written in Schwarzschild coordinates as ^fgdx 1 . . . dx 4 = r 2 sm(0)drd9d(f)dt. 

The estimate in (|1.5p reveals the importance of the photon sphere, r = 3M, in this problem. The null 
geodesies which orbit at r = 3M appear to have some similarity to closed geodesies on a Ricmannian manifold, 
since their spatial coordinates are periodic in timcE However, because of the curvature of both space and 
time, it is not clear how closely this analogy holds. Physicists are well aware of various strange effects caused 
by the photon sphere, for instance, an infinite number of visible images are created by gravitational lensing 
around the photon sphere |15j . 

The theory of dispersive waves on manifolds is still developing, but dispersion is known to be weaker in 
the presence of closed geodesies. Although the exponents differ, in K n+1 , the wave and Schrodingcr equations 
have similar L 1 — > L°° pointwise in time decay rates and similar Strichartz estimates |51 130j. In the absence 
of closed geodesies, global Strichartz estimates have been proven for both the wave [35J and Schrodingcr 
[21j . For the wave equation, because of finite speed of propagation, the local in time theory is the same as in 
jjn+i. nowev er, for the Schrodingcr equation on a compact manifold with closed geodesies, on a finite time 
interval, Strichartz estimates hold, but with an e loss of regularity [7]. Similarly, while an estimate like (|1.6p 
holds for the Schrodingcr equation in R" +1 (where it is known as a local smoothing estimate, since control 
on the left is by an H 1 ! 2 norm for the Schrodingcr equation), there must be an additional e derivative loss 

1 Morawetz originally used a somewhat different estimate, which was modified to prove local decay estimates of the type we 
use in R n 1231 and later for radial data on the Schwarzschild solution 1221 . Morawetz also introduced what we refer to as the 
conformal estimate in R* 1 , which some authors refer to as the Morawetz estimate. 

2 We do not claim that the behavior of null geodesies tangent to the photon sphere is analogous to closed null geodesies on 
a Lorentz manifold, which violate causality and which are not present in the Schwarzschild manifold. 
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in the presence of closed geodesies [T3]- For the wave equation in the presence of closed geodesies, it is not 
possible for an estimate like (|1.6[) to hold if it controls derivatives in all directions, since, through a geometric 
optics approximation, wave equation solutions in the energy space can be shown to stay arbitrarily close 
to the path of geodesies [53]. Since we are working with commutator methods, which typically can also be 
applied to the Schrodinger equation, we do not expect to be able to avoid an e regularity loss in the presence 
of closed geodesies. 

The analogy we wish to consider is to a wave equation on a three-dimensional, warped product manifold, 
Wl = R x 2U, with metric ds 2 = dr 2 + r(r*) 2 dtu 2 and 2U compact. We make several assumptions on r(r*), 
explained in section [21 but essentially that r has a unique minimum, which corresponds to a single closed 
geodesic surface, and that r grows at least linearly but at most polynomially. This Riemannian metric defines 
a Laplace- Beltrami operator, which we can use to define a wave equation, 

-d 2 4> = -A m fa 0(0) =4>o, S t 0(O)=0i. (1.8) 

Introducing <p = r<f) again, we get an equation of the form (| 1 .3ft . The analysis we develop for the wave 
equation on the Schwarzschild solution applies equally well to this problem both strengthening the analogy 
between the photon sphere and a closed geodesic on a Riemannian manifold and proving new estimates 
for solutions to wave equations in the presence of closed geodesic^. We are able to prove the analog of 
proposition ll.il In fact, for the defocusing, semi-linear wave equation with exponent 8/3 < p < 3, 

-0 2 ^ = ~ A^ + IC" 1 ^ 0(0) =0o, 0t0(O)=0i. (1.9) 

our space-time integral estimates are sufficiently strong to close a boot strap argument and extend our result 
to the small-data, semi-linear case. 

Proposition 1.2. If <f> is a solution to the linear wave equation (jl.8p on a warped product manifold, s > 0, 
and conditions [7] to given in section \^ hold, then 

\4>\ 4 d 3 voldt < OE c [0o,0i] + £[£ £ 0o,£ £ 0i]) 

oo JWl 

provided that the norms on the left are sufficiently small. The operator L denotes the operator square of 
(1 — Agu). For 8/3 < p < 3 and sufficiently small initial data, in the norms appearing on the right hand 
side, the same estimate is true if (p is a solution to the semi-linear wave equation (|1.9p . 

Because this result replies upon commutator methods, it can also be extended to the Schrodinger equation 
on a warped product manifold, using the methods in this paper and in [5] and |22j . We will deal with this 
elsewhere. 

We begin in section [2] by listing the conditions on the terms in the wave equation (jl.3[) which allow us 
to analyze it and explaining how to transform the wave equation on the Schwarzschild solution or a warp 
product manifold to fit this description. In section [31 we define the energy and conformal charge. In section 
21 we prove a Sobolev estimate, an estimate on weighted L q norms in terms of the conformal charge, and 
a Hardy type estimate. Following this, in section O we introduce an analog of the Heisenberg relation for 
the wave equation, which we use in section [6J to prove the local decay estimate (|1.6|) . In section [3 we define 
and outline the use of the phase space observable, T, which allows us to improve the local decay estimate 
by a factor of L 1 " 5 . In section [51 we then localize this inside the light-cone to get additional decay in time, 
which allows us to close the estimates and prove Propositions 11.11 and 11.21 In the appendix, we present a 
brief version of the arguments in [5] which allows us to estimate the partial phase space observables, T n ^ m , 
which make up T. 



2 Setup 

In this section, we state the form of a wave equation to which the rest of our analysis will apply, and we 
show that the wave equation on both a warp product manifold satisfying certain conditions and on the 

3 As noted in [6], the higher regularity analysis there also applies to warped product manifolds with a single closed geodesies 
under certain conditions; although, the necessary conditions were not given. 
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Schwarzschild solution can be transformed to this form. In addition to the differential equation, we will also 
introduce function spaces. 

We will consider a wave equation of the form 

-d?0 = -d?J + V<P + V L (-A w )<j ) , 0(1) =0o, 0(1) =0i- (2.1) 

for cj)(t, r^to) with (t,r*,ui) G R x DJl = M x R x 2U and 2U a two dimensional, compact, smooth manifold. On 
2B, the Laplace- Beltrami operator, (— Agu), has a spectrum, a{— Agu), which is discrete [2]. Following the 
standard terminology on the sphere, we refer to the eigenspaces as harmonics. We index the eigenvalues in 
increasing order by I, and denote the projection onto the corresponding cigenspace by P/ and the eigenvalue 
by I 2 . Thus, on the sphere, we have I 2 = 1(1 + 1). We will also consider the nonlinear wave equation, 

-9 t 2 = -a r 2 > + y0 + v/ L (-A 2D )0 + /(r,)|0rV, 0(1) =00, 0(1) =01 ■ (2.2) 

We denote the nonlinearity F'(|0| 2 )0, where 

I0| P+1 



p+l 



We treat the linear wave equation (|2.ip as a special case of the nonlinear wave equation (|2.2p with / = 0. 

Since these equations are extremely general, we only consider the case where the following assumptions 
hold: 

1. The potentials, V and Vl, are smooth functions and depend only on the radial coordinate, r*. 

2. The potentials are non-negative: V > 0, Vl > 0. 

3. For each I 2 in the spectrum of (— Agjj), for the effective potential Vi = V + 1 2 Vl, 

(a) Vi has a unique critical point, which is a maximum and occurs at r» = (aj)*, 

(b) Vi vanishes linearly at r„ = («;)*, 

(c) the potentials V and Vl also satisfy conditions I3al and I3bl with (ao)* and (aoo)* respectively. 

(d) the (a;)* converge to (a M )». 

4. The trapping terms, IV + r*V and 2Vl + r t V£, are positive only in a compact region. 

5. The weight V[/Vl is bounded above by a constant, and the angular potential decays at least like r~ 2 : 
V'JVl <C,V l < Ca + r 2 )" 1 . 

For non- linear problems, we require the following additional assumptions: 

6. The weight V^/Vl decays like r" 1 , and the angular potential decays at least like r~ 2 : V[/Vl < 
C(l + \n\)-\ V L <C{l + rl)-\ 

7. The non-linear coefficient is a smooth function which depends only on the radial coordinate and is 
nonncgativc: / > 0. 

8. The non- linear trapping term, 2/ + ?'*/', is positive only in a compact region on which / is strictly 
positive. 

9. The quantity fV^~ p ^ 2 is bounded. 

Note that condition [S] follows immediately from condition [51 Note also that the existence of a unique critical 
point which is a maximum, from condition I3al guarantees that the potentials are bounded. Furthermore, 
from condition 3] on the trapping term, the quantity r 2 V and t 2 Vl are bounded, so the potentials must 
vanish rapidly as r* — > ±oo. For the local decay, angular modulation, and phase space analysis techniques 
used in sections [6] and [71 conditions [Tl - [3dl are sufficient. 

We will work with the measure d 3 /x = dr*dui. Unless otherwise specified, L q norms are with respect to 
this measure, \\?p\\ refers to the L 2 norm, and inner products are taken in L 2 . There is a Fourier transform 
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on £DT = R X 2B denned by decomposing in the eigenbasis for (— A%j) and then applying the one dimensional 
Fourier transform on each eigenspace. The Schwartz class, §, is defined to be the set of pairs (V'O) V'l) f° r 
which 

Vi, j, k G N,a G {0, 1} : iC^,* : V(r*w) G SOT : |r^, V^aCr*. w)| < C Q , 4J , fe . (S) 

Abusing notation, we will use ^(£) G § to mean (il>(i),ij>(t)) G S. In general, we will use <p to denote a 
solution to (|2.2p . V' or ?/>(i) to denote a function in § in the abused notation, and (v,w) to denote a pair in 
S. 

In our analysis, the function <p will be related to a function, <p, which is the solution to a geometric 
equation. These are related by <p = r(r^)<p for some weight function r(r„). 

We now consider the wave equation on a warp product manifold. We take the manifold to be (r*,o;) G 
dJl = R x 2B, where 2U is a two dimensional, compact, smooth, Riemannian manifold with the metric du 2 . 
We take the metric on 9Jt to be 

ds 2 =dr\ + r(r*) 2 dw 2 . 

This defines a Laplace- Beltrami operator on 9JI, and we consider the nonlinear wave equation for (f>, 

-d 2 t j> = -A m j> + \^- 1 ]>. 
Under the transformation <fi = r<p, this becomes 

-d 2 cf> = - d 2 r _^ + Vcjy + V L (-A W )<P + /F'(|0| 2 )0 

with 

V =r"r -1 , 
V L =r~ 2 , 
f =r X - p , 

p+l 

This is of the form (|2.2|) . If r and r" arc positive then conditionsHJEUll and[9]are clearly satisfied. Conditions 
El IH[6l and[8]arc, therefore, the only restrictions. From the example r = (1 + r 2 ), it is clear these conditions 
permit some manifolds with closed geodesies. If we consider the restriction on Vl to be the central ones, then 
condition [3] says there there is a unique, closed geodesic surface which is unstable. If r behaves polynomially, 
with r = rl 1 + 0(r^ 1_1 ) and r' = r, 1 + (^(r* 1 ), then conditions [4] and [6] follow if r grows faster than 
linearly, pi > 1 (and will also hold if pi = 1 provided the lower order terms give condition HJ. This is also 
sufficient for [8] to hold. We work with the measure dr*d 2 Lo for <fi and with r 2 dr*d 2 uj for <j>. This gives an 
isometry of L 2 spaces. 

We now turn to the wave equation on the Schwarzschild manifold. Since there is curvature in time as 
well as space, the wave equation does not arise from a warp product of the form previously considered. 
The exterior region of the Schwarzschild solution is given by (t, r, uj) G R x (2M, oo) x S 2 with the metric 

ds 2 = - (1 - 2M/r)dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 duj 2 . 

We introduce the Regge-Wheller tortoise coordinate, r*, defined by 

dr 

— =1 - 2M/r, r(0) =3M. 

dr» 

It is sometimes useful to work with the explicit expansion, 

(r — 2Af\ 

U = r + 2M log I — J - 3M + 2M log 2. 
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This transforms the exterior region to (t, r*, lu) e R x Wl = K x K x S 2 . 

The Lorentzian metric defines a d'Alembertian operator and thus the wave equation (|1.2[) for </>. Using 
the coordinate and the transformation cf> = r<f>, the wave equation takes the form (|2.ip with 

2M 

V=—(l-2M/r), 
V L =l(l-2M/r), 

Once again, it is clear that conditions [U [2l and [5] arc satisfied, since we are working in the region r > 2M. 
We show conditions [3] and [4] hold below. We use the measure dr*d 2 LO for <fi and use the natural measure 
r 2 drd 2 iodt = y / —gdrd 2 ojdt = (1 — 2M / r)r 2 dr*d 2 u>dt when integrating (j) in space-time. 

If we were to introduce a semi-linearity into the wave equation, this would translate into a term fF' ( \(f>\ 2 )(j> = 
(1 — 2AI/r)r 1 ~ p \(j>\ p ~ 1 (f> in equation (|2.2[) . This clearly satisfies conditions [7] and O For p > 3, it is not hard 
to show [5J Unfortunately, condition [5] does not hold towards the event horizon, r — > 2M, which prevents 
us from getting a good bound on ||r -1 (/)|| in the Sobolev estimate and from controlling the growth of higher 
Sobolcv norms. 

We first show that on each spherical harmonic, the effective potential has a unique maximum. In [5], 
we showed the same holds true on the more general Reissner-Nordstr0m manifold, but, for simplicity, we 
restrict to the Schwarzschild manifold here. 

Lemma 2.1. Fori > and M > 0, the effective potential Vi has a unique critical point, which is a maximum. 
The maxima, (a;)*, converge to 0. 

Proof. From the effective potential, we can compute its derivative, which has two roots, ai_±. 

Vi =V + l 2 V L , 

V{ = - 2r- 5 (l 2 r 2 - (P - l)3Mr - 8M 2 )(1 - 2M/r), (2.3) 



(I 2 - 1)3M ± J {I 2 - 1) 2 9M 2 + 32l 2 M 2 
«M= ^ ■ 

For I = 0, the roots are and 8M/3. For I > 0, ai t - is negative, and at = is larger than 2M. At r — ai, 
V{ goes from negative to positive. From (|2.3[) . it is clear that ai — > 3M. The corresponding values of r* are 
denoted (a/)* and converge to (ctoo)* = 0. □ 

By direct computation, we now show that the trapping terms 2U + r*U' and 2Vl + r*V' L are positive 
only in a compact region. 

2AL2M , 2M.2M, 
2V + nV =2(1 - — )— + r*(l - — )_(_3r + 8M) 

= 2M 2M a _ _ 
r° r 

2 2Af 1 
2V L + r t V[ =-^(1 - — ) + n-j(-2r + 6M) 

=4(l_f^)(r»-r.(r-3JW)) 

Since for large r, r* ~ r + logr, the terms 2r 2 — r*(3r — 8M) ~ — r 2 and r 2 — r*(r + 3M) ~ — rlogr 
both diverge to negative infinite. For r-„ going to negative infinite, 3r — 8M < and r — 3AI < 0, so that 
— r»(3r — SM) and — r*(r — 3Af) both diverge to negative infinite like a constant times r*. This proves that 
the trapping terms 2U + r*V' and 2Vl + f*V£ are positive only in a compact region. This shows condition 
HI holds. 

From condition[Sl Schwartz class functions and all their derivatives must vanish at the bifurcation sphere. 
The normalized radial derivative is (1 — 2M/r) 1 / 2 d rt ip. Since r — 2M and r* are related exponentially as 
r — > 2A/, the polynomial decay rate for r* derivatives does not imply decay, or even boundedncss, of r 
derivatives or normalized radial derivatives, which more accurately describe the behavior of ip as r — > 2M . 
Thus, functions like (r — 2M) e satisfy condition [5] as r* — > — oo (some other factor which decays as r — > oo 
can then be applied to enforce the condition in the other direction). 



8 



3 Methods for pointwise decay 



In this section, we summarize the use of various densities, in particular the energy density, which we have 
previously used to study the wave equation on the Schwarzschild space [5]. This is an application of very 
commonly used techniques. 

The energy, radial momentum, and (T2U valued) angular momentum densities are defined by 

e[v, w] =^(w 2 + v' 2 + Vv 2 + V L X7 w v ■ X7 w v + /(r»)F(u 2 )), 
p*[v, w] =wv', 

If (j) is a solution to the wave equation, these satisfy 

=d t e[<j),(j)] - d r „p* [</>,</>] - Vac • p w [</>, <j>], (3.1a) 
=d tP *[<t>, 0] - \d r , [# - <t>'4>' ~ <f>V<f> - Vau^. • V L Vw^ - /(r*)F(0 2 )] (3.1b) 

- Van ■ (4>'V L Vw^>) 

- \<t>V'<t> - \^ W <P • Vlymt - \f{n)F(\<}>\ 2 ). 
From integrating equation (|3.1a[) , we have that the total energy, 

E[<f>,4>] = f e[<t>,j>]cPn= f \{^ + ^ 2 + V4> 2 + V L V w cp-V w ^ + f(r r )F(cf> 2 ))d^ 

is conserved. Since the energy density dominates the 1? norm of the time and radial derivatives of <f>, we 
can use many of the basic arguments which were developed in K n+1 ; however, since the weight, Vl, may 
have more complicated dependency, we will only be able to use the angular derivatives through a Sobolcv 
estimate and the control in condition or [SJ 

We now define the conformal charge density in formal analogy to the corresponding density in K n+1 , 
where it is generated by a symmetry of the wave equation. In the present situation, it does not come from 
a symmetry and does not lead to a conserved quantity. We define the conformal charge to be the integral 
of the conformal charge density plus the energy density. This slightly strange definition is taken to simplify 
the statement of certain estimates, for example estimate (|4.3j) . 



e c [v, w] =(t 2 + r 2 )e[v, w] + 2tr*p*{v, w] 

Ec[v,w] = / ec[v,w] + e[v, w]d 3 /i 
Jm 

It has long been known that the conformal charge density can be written as an inherently positive quantity 

ED- 

e c [v,w) =^{t-n) 2 (w-v') 2 + j(t + nf(w + v') 2 (3.2) 
+ \{< 2 + rl)Vv 2 + l -{i 2 + r 2 )V L (V w v ■ V w v) + \(t 2 + rl)f(r*)F(v 2 ). 

In fact, the control on t 2 times the angular derivatives is the central element which allows us to prove 
decay of LP norms through the Sobolev estimate. The conformal charge can be shown to control the L 2 
norm, and, interpolating this with the Sobolcv estimate provides the following estimate, which is proven in 
lemma B~2l for t/j 6 S, 

f(r*fVjT V| 3(p " 1) d 3 /' KCEf^t-P+i. 

m 
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We now turn to estimating the growth of the conformal charge. Introducing further the Lagrangian and 
the shear densities 

l[v,w] = ^(w 2 - v' 2 - Vv 2 - V L W w v ■ W w v - f(r*)F(v 2 )), 
T**[v, w] =v v , 
T* u [v, w] =u'VjrV2n«, 

the relations (|3.1ap - (|3.1b|) can be rewritten as 

=d t e - d rt p* - Vaj • p u , 

Unless otherwise indicated, we will evaluate all densities on a solution and its t derivative. (<fi, <f>). From this 
presentation, and integration by parts, it is easy to show that 

d t ee =2te + 2r*p* 

+ (t 2 +rt){d r ,p* + Vau -Pu,) 

+ 2tn(d r J + d r ,T„ + Vau • T w + ^V'tf + ^Vl\V w <j)\ 2 + ^/'(n)F(0 2 )), 

-E c = / 2t(e - I - T 1 **) + tr*{V'<t> 2 + V[\V w 0\ 2 + f'F{cf> 2 ))d 3 fi 
at J m 

= [ t(2V + r*V')<j } 2 +t(2V L +r»Vl)\V w ^\ 2 +t{2f + rJ')F(<j) 2 )d 3 fi (3.3) 
Jot 

To control the growth of the conformal charge, it is sufficient to control a space-time integral which 
is localized. Introducing Xa to denote a smooth, positive, compactly supported function which dominates 
2V + r*V, and 2Vl + r*V[, for linear solutions, 

E c [<t>A}{t) <C [ r||xa^|| 2 dr + £:c[<M(0) (3.4) 
Jo 

For nonlinear solutions, the integral J Q r j m xF(\cj)\ 2 )d 3 fidr must also be controlled where \ also dominates 
the nonlinear trapping term, 2/ + r*/'. 

4 Three weighted norms controlled by derivatives 

We present three lemmas in this section. The first is a Sobolev estimate on DJl, representing either a warp 
product manifold of the type we consider or a spatial slice of the exterior region of the Schwarzschild solution. 
The second is an L p estimate in terms of the energy and conformal charge, which we need in estimating the 
nonlinear contributions to the growth of the conformal energy and the growth of the higher angular Sobolev 
norms. The third is a Hardy or Poincaree estimate, which controls a weighted L 2 norm by the sum of the 
L 2 norm of the radial derivative and the expectation of the potential. 

We now show that a Sobolev estimate holds both for warp product manifolds of the type we've considered 
and for the Schwarzschild solution using condition The main feature of this estimate is that it separates 
the role of the angular derivatives from that of the radial derivatives. This is important, because we will 
later show that the angular derivatives decay 

Lemma 4.1. For ip € S, 

f VM e d 3 n<C([ \d r ,ij\ 2 + \(Vl/V L W\ 2 d 3 v)([ ^(ivrf + H 2 )^) 2 - (4.1) 
Jm Jm Jm 
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Proof. We reproduce the H 1 <^-> L 6 estimate presented in [5] and following the standard methods for R 3 |16j . 
We begin by proving a W 1 ' 1 ^ L2 estimate and define 

= / Idr.V'oO"*) w)|dr* 

JM 

We use f/'o to temporarily denote a function with unspecified relation to ip. From the Fundamental Theorem 
of Calculus, the Cauchy- Schwartz estimate, and the W 1 ' 1 (W) L 2 (2U) Sobolev estimate, for ip G §, 

|^o(^*,w)|5 <7i(w)l?|V>o(r-*,w)| 

|Vo(r*,w)|5dw <( / h(uj)duj)?{ \^ (r*,uj)\ 2 duj)^ 
s 2 Js 2 is 2 



<(/ |9r,^o|d>) 2 / iVauVol + 
ion Js 2 

Integrating in r* against the unspecified weight h 3 gives 

h^\M i d 3 ^<{[ l^.^old 3 /*)' / ^(IVauVol + IVoDrfV 
931 Jot Jot 

We substitute ^0 = ^l^l 4 an d apply the Leibniz rule and Cauchy-Schwartz, 

h 2 \i/j\ 6 d 3 n <4(/ h 2 \^\ 6 d 3 ^(f \d r ,ip\ 2 + \h'h- 1 iP\ 2 d 3 fJ l )i( f h\V w ^\ 2 + h\ifd 3 ^. 
ot Jot Jot Jot 

We now take h — Vl- □ 

The following estimate shows that certain L p norms are controlled by the energy and conformal charge 
with decaying powers of t. The decaying powers of t provide the decay estimates necessary to show that the 
nonlinear contributions to the growth of the conformal charge and the higher angular Sobolev norms are 
small. 

Lemma 4.2. For ip S §, under condition^ 

VM 6 d 3 fi<CEc[>P,^} 3 t- A (4.2) 

OT 

Further, under condition^ 

[ fV L 1 \^- l U 3 i i<CE^ 11 t 3 ^ + i\ (4.3) 
Jot 

Proof. We begin by explaining how to control the inhomogeneous terms, involving no derivatives, on the 
right of equation (|4.1|) . This involves controlling the L 2 norm by the conformal energy, following [18] , and is 
similar to the Poincare estimate (|4.7[) . 

Focussing on the time and radial terms in the conformal charge, 

2E c [4>,i>] >\\\{t-r){^- d r M\ 2 + \\\(t + r*)(V3 + d r ,4>)\\ 2 = N + r^M 2 + \\V4 + td r ,^\\ 2 . 
Adding and subtracting — hr 2 (l + r 2 )~ 1 -0 and hr^t{\ + r 2 )^ 1 ^, and expanding, we have 

E c >\\\4 + r*d r ,Tl> - hrl{\ + r 2 )~Vli 2 + + td r ,^ + hr,t(l + r 2 )"Vl| 2 

- (d r ^J, (t 2 - r 2 )r,(l + r 2 )- 1 /^) - \{h^, (t 2 + r 2 )r 2 (l + r 2 )" 2 /^). 



Integrating by parts and noting that r* is differentiable, the last two terms can be written as 

2 {l + riy i + n 
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For h E (—1/4,0) and |r*| > 2, this weight is strictly positive, and using the potential V to control small r*, 
we have 



EcbP, i>] >\\\{t - r){i> - d r .M\ 2 + \\\(t + r.)(V- + ^.VOf 



+ \ I (t 2 +rl)ViP 2 d 3 [i 
t 2 + r 2 



j m i -t- r„ 

>CU\\ 2 - (4.4) 

By conditional wc can estimate the (V[/Vl)iP 2 term by ip 2 and the Vl^ 2 term by (1 + r 2 )~ 1 xp 2 . Since 
E < E c , 

[ Vl\i>\ e d 3 n<CElt-\ (4.5) 
Interpolating this against (|4.4|) . for 2 < a < 6, 

V^ip a d 3 fji <CE^r a+2 . (4.6) 

93! 

3(p-l)-2 

To conclude, we consider the case a = 3(p — 1). From condition [9l 2 > f V L , and the desired 

result holds. □ 



Remark 4.3. T/ie estimate (|4.4p requires that ip — > as r* — > ±oo. These are natural boundary conditions 
on a warp product, but not entirely natural as r* — > — oo on i/ie Schwarzschild solution. The Sobolev estimate 
in lemma \jn\ still holds for functions which do not decay as r* — > — oo, since t/ie radial integration in the proof 
can be taken from r* — > oo, where there is decay. The estimate (|4.4[) does require decay at —oo; however, by 
applying the same argument to the region r* > 0, and using the contribution from the potential term in the 
conformal energy, (|4.4p can be replaced by J m (l - 2M/r) 2 ip 2 d 3 fj. < E c . This can be interpolated against the 
Sobolev estimate to control f m (l - 2M/r) 2 r" ff+2 f d 3 /*- With cr = 3(p - 1) and f = (1 - 2M/r)r 1 - p , th 



e 



integrand is exactly f 3 V L 1 ip 3 ^- p 1 >. Thus, this lemma would hold for the semi-linear wave equation on the 
Schwarzschild solution. It is condition^ not\Q which prevents the main result from holding. 



Remark 4.4. Under condition^ it would be possible to estimate the term \\r ip\\ by the Hardy inequality 



below, so that the weighted L a norm would be controlled by CE a ^E a J' A t CT+2 



The following Hardy or Poincare estimate controls a weighted L 2 norm by the L 2 norm of the radial 
derivative plus the expectation value of the potential. We use the contribution from the potential, since our 
method is to reduce something resembling a one dimensional wave equation, and the radial derivative alone 
is not sufficient in one dimension. The purpose of this estimate is to show that 7 is energy bounded. 

Lemma 4.5. For ip € 8, 

\\(l + rlyh\\<CE[ij^}. (4.7) 

Proof. This argument is taken from [5J. We start by considering a one dimensional, bounded, H 1 function, 
which we will also denote by ip but with argument x, and then integrate the derivative of (1 + |x|) _1 '0 2 
against the weight sgn|a;| on a bounded interval. 

-W>( Cl ) 2 + V(-ci) 2 )-2V;(0) 2 = / 1 S gn(x)d x (-^- 1 ^(x) 2 )dx 



1+ ci ' rv ' ' rv ' y_ ci ° v ' v i+ x 



ci 
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Isolating the integral of the weighted term, we have 

;dx = - —L-^( Cl ) 2 + V>(-ci) 2 ) + 2^(0)^ + / sgn(x)— ^^d^dx 



- Cl (l + M) 2 l + | Cl | v ^^ ^ rw J_ Cl aK 'l + \x 

<2</>(0) + (J 1 (1+ 2 | X | )2 V' 2 ^) I ( 1 1 {d^fdx) 1 -. 

Since the weighted norm appears quadratically on the left, but only linearly on the right, taking the limit 
as ci — ► oo, 

{l + \x\y dx ^ C (^(0) 2 + I9^| 2 ^)- (4.8) 

On 9Jt, we now generate a family of one dimensional functions by fixing each uj and t variable and by 
translating in the radial direction by a quantity co, that is taking x = ro — cq. Integrating ()4.8|) over a 
compact interval of cq values, 

^* )2 r dr, < C( [°° V^ 2 (r*) + |^V| 2 dn). 



(i + N) 2 



Integration in the angular variable, and estimating the right hand side by the energy gives the desired 
result. □ 



5 Heisenberg relation 

We now continue with an approach based on the interpretation of the Schrodinger equation as describing 
the evolution of a quantum system. This allows us to think of certain operators as corresponding to physical 
observables. We begin by writing the wave equation in terms of an operator, 

4> + H(j) = 

where H is the Hamiltonian defined in terms of the following pieces 

H =Hi + H2 + i?3, 

H\ = — d 2 r , 

H 3 =V L (-A W ). 

The Hamiltonian, H, and it pieces, Hi, are self-adjoint. From this, we have a Heisenberg relation describing 
the time evolution of a skew expectation value of an observable. 

Theorem 5.1. If A(t) is a differentiable family of self- adjoint operators with dense common domain, and 
if (p G S is a solution to the non-linear wave equation (|2.2p . then 

j^M) - (<M$) =fa> [H,A]<f>) + (<t>,A<j>) - {^M) + (0,A/(r»)F'(0 2 )0) - (/(r»)f'(</) 2 )<M</>), (5.1) 

where (0, [H, A]<fi) is understood to be the quadratic form defined on D(H) D D(A) by (Hcj>, A<f>) — (Acj>, H<f>) . 
This result follows by direct computation. 

Our method for proving local decay and the phase space estimate is to find a propagation observable. A 
simple example of a propagation observable which majorates an operator G is a time independent operator 
A for which 

[H,A] = G*G. 
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If A is energy bounded, in the sense that ||^4</>|| 2 < E[(j),<j)\, then from integrating the Heisenberg relation 
(|5.1[) . we have a space-time integral estimate on G: 

J \\G<f>\\ 2 dt< J j t ((<f>,Acj>)-(t,Act>))dt<CU\\\\Act>\\ <C£7[0,0]. 

The notion of propagation observable is slightly broader in that it allows various lower order terms. 

Definition 5.2. Given a pair of possibly time dependent operators, A andG, with D(A)nD(G)nD(H) dense 
in L 2 (9JV), the operator A is a propagation observable which majorates G if, for all ip G D(A)nD(G)l~)D(H), 

(V, [H, A]ip) + Aj>) - (V>, M) >(Gip, &P) + (V, Ril>), 
where R is a sum of lower order terms (with respect to G), Ri, for which either 

1. Rj is a bounded operator and if <p G S is a solution to the non-linear wave equation (|2.2[) . 

J (<p,R l cp)dt < C 

2. or Ri is an operator with D(Ri) D D(G*G) dense and for all ip € D(G*G), 

\{^,Ri J )\<C^,(G*G) 1 -^) 

We make the following additional remarks. The notion of lower order depends on the operators A and 
G. HA is an anti-self-adjoint operator, then the same Heisenberg relation (|5.ip holds. If H and A take § to 
§ then the commutator acting on S is the operator given by 

[H, A] = HA — AH. 

When making estimates involving operators B and R, with S dense in their common domain, we will 
abuse notation by writing 

B > R 

to mean that for ip € S 

0>,Bip) > {^,Ri>). 

6 Local Decay 

The goal of this section is to prove bounds on weighted L 2 space-time norms of solutions. We begin by intro- 
ducing a propagation observable, 7, show it majorates a weighted quantity, and conclude with a Gronwall's 
type argument to integrate the Heiscnberg-type relation. These estimates are proven using an eigenfunction 
decomposition analogous to a spherical harmonic decomposition, we refer to this as a harmonic decompo- 
sition. We show these estimates have a uniform nature in the harmonic parameter, so that we can recover 
an estimate for general u. In later arguments, the contribution from H2 is controlled by the local decay 
estimate, and we can use a uniform multiplier. 

In K 3+1 , the radial derivative can be used as a propagation observable to control the time integral 
of |it(i,0)| 2 . This can be thought of as a weighted space-time integral with the 8 function as a weight. 
Essentially, our propagation observable is a smooth version of the radial derivative, which leads to a smooth 
weight. 

This follows [5^1 [3]. Condition [3J on the effective potential, Vi, is the main condition used in this section. 
We recall the effective potential on each spherical harmonic, given by 

Vi =H 2 +H 3 = V + l 2 V L . 
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If necessary, we will use P; to denote the projection on to the harmonic with parameter /, but typically, in 
this section, we will assume that u has support on only one harmonic. In this case, <p satisfies 



(/)-(/)" + Vi4> = o. 

From condition I5al Vi has a single critical point, which is a maximum and which we denote (a;)*. 

The Morawetz-type operator is defined on each harmonic in terms of the radial coordinate r*. The 
definition involves a weight, gi M , which is defined for a > 1 so that </; jCr remains bounded. There is also 
a dilation factor b which we will later choose sufficiently large, independently of Z, so that 7 majorates 
(1 + ft) • Wc will prove this on each harmonic, and then show there is a uniform lower bound on these 
estimates, to get an estimate uniform in I. 

Definition 6.1. Given a > 1, b > 0, and I € N, the Morawetz-type multiplier 7; jCr is defined by 

Mr, 1 

9lAn) =Jo WTW dT 

=9iA r *) d r, + ^9'i,*(r*) 

In all cases a will not vary so the notation 7 = 7; >cr and g — g\ a will be used. When working on more than 
one harmonic, we also use 7 to denote 

l 

We note that as / varies, the observable 7; j(T is simply translated. Since the (a;)* converge, the weights 
gi t a has a limit in L°°. 

We now show that the commutators uniformly bound a polynomially decaying term, for an appropriate 
choice of b. This value of b will always be used in 7j j(T from now on. There are additional positive terms 
which are also dominated, but we are better able to take advantage of them in the later sections using a 
slightly modified multiplier. 

Lemma 6.2. If a > 1 and M > 0, then there is a choice of b for which there is a constant C such that for 
all ^gS, 

(ip,[H,l]Tp) >C(V,(l + r2)-f-V) + C(9 r ,V,(l + ^)-t9 r ,^). (6.1) 

Proof. We work on a single spherical harmonic and sum at the end of the argument. The commutator can 
be computed exactly as 

\H,l] = ~ \<&.gd r . + d 3 rt g - gd^ - d r ,gd 2 r J - (gd^Vt - V l9 d r ,) 
= -2d r ,g'd r , -\g"'-gV{. 



We will use x to denote r* — (a;)» 



2' 



9 ' = (TT%[r (6 - 2) 

_ -6 2 <7 S gn(a:) 

g -(i + b\ x \y+i (b - 6) 



To estimate the contribution from the S function, we introduce a smooth, weakly decreasing, compactly 
support function, \ which is identically one in a neighborhood of x = (this is unrelated to any other cut off 
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function used elsewhere in this paper). Using integration by parts, and temporarily treating ip as a function 
of x, 



0=1 d Tt (xxtp )dx 

X?p 2 dx + J x\ip 2 dx + J 2xx>pd Ts> ipdx, 
Xi> 2 dx < j (\x X '\ + x 2 X )tp 2 dx + J x(d r ,ip) 2 dx, (6.5) 
^(0) 2 =- / d r ,(x^ 2 )dx 

poo 

x'ip 2 dx- / 2x4>d r ,ipdx (6.6) 
Jo 

poo />oo 

</ W\ip 2 dx+ Xl p 2 dx+ x{d r ,i>) 2 dx. 



Applying a symmetric argument for (— oo,0] and substituting the result for J xip 2 dx, we have 

2^(0) 2 < J (Ix'l + \x X '\ + x 2 xWdx + J 2x{d r Mdx. (6.7) 
We use a similar integration by parts argument to estimate (1 + b\x\)~ a ~ 2 il}. 

= j dr, [g"^ 2 )dx 
= I g"'ip 2 dx+ / 2g"ipd r ,i^dx 



1} T *jj 2 dx <2ab 2 yj(0f 



(1 +b\x\)° 



+2 



,3 <r(a + l) ,,, \ * ( f 2 * , 2 (dr^) 2 



67 — ,, -x ^ dx / b , v dx (6.8) 

(l + b\x\Y+ 2 J \J a + l (i + b\ x \y 



Substituting equation 
a(a + l) 



(l + b\x\) 



er+2 



iP 2 dx <ab 2 / (Ix'l + \x X '\ + x 2 X )Vdx + ab 2 \ 2 X (d r ,tp) 2 dx 



1 f. a(a+l) - f 2a , 2 (d r ^ 2 , 

2 7 (1 + 6 a CT+2 7 <r + l {l + b\x\Y 



l) ,i, 2 dx <2ab 2 f (Ix'l + \x X '\ + x 2 xWdx + 2ab f (— \ rT - + 2b X )(d r ^) 2 dx 



(l + b\x\Y+ 2 ~ 7 7 cr + 1(1 + 6|x|)° 

(6.9) 

The commutator term can now be estimated 
2g'(d r ^) 2 -±g"'iJj 2 ~gV l 'iP 2 dx 

{d r ,^) 2 dx + ab 2 yj(Q) 2 - \ ( ^ a{ °\ l l^ 2 dx - gV{^ 2 dx 



{l + b\x\Y rw 2 7 (l + 6|a;|) CT + 2 

r 2b 

- J ^TbWY {drJj)2dx + ab2m2 

- °b 2 J (Ix'l + \x X '\ + x 2 xWdx -abj ( - | - - 1 + 2b X )(d r ^) 2 dx - J gVtfdx 

>2b f ( T-TTT-rrr iw+ 2 - ^ x )(9 r ,V) 2 ^ + / (-^ 2 (lx'l + \*x'\ + * 2 x) - 9Vl)^ 2 dx + ab 2 ^(0f 



cr + l {l + b\x\Y+ 2 



(6.10) 
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The term abtp{0) 2 is positive and will be ignored. The term ((7 + l)~ 1 (l + fo|a:|)~' T ~ 2 —abx is independent of I, 
so we can choose b sufficiently small so that this is bounded below by c(l + b\x\)~ a ~ 2 . Since V{ — V + V[l , 
and the (a;)* converge, the value of V/' is uniformly bounded from below at (a/)*, and b can be chosen 
sufficiently small so that the quadratically vanishing quantity — er& 2 (|x'| + \ x x'\ + x2 X) ~ 9^{) is bounded 
below by c6 2 (|x'| + \x\'\ + x 2 \). l uniformly in Thus there is a uniform constants, Ci, for which, 

>ci J {1 + l\ x \ r ( d r^) 2 + ° b2 (\x'\ + M + x 2 xWdx 
Applying the estimate (|6.9p . there is a uniform constant, C2, such that 

2g'(d r ^) 2 -^g"'^ 2 -gV l ^ 2 dx 

^ I (i + b\x\r + ^ dx + C2 / (TtW (9 ^ )2 ^ (6 ' n) 

Since the converge, there is a uniform equivalence between (1 + 6|a:|) _<:r_2 and (1 + r 2 )~( cr/ ' 2 ' 1 . Inte- 

grating over the angular variables gives 

(i + r 2 )2+ i {i + n)2 

□ 

For the local decay result, equation (|6.1|) is sufficient, but in the phase space analysis, wc will need more 
contributions from —gV[. We now prove the local decay estimate. Since 7 is self-adjoint and preserves the 
class of real valued functions 

(^7V)-(V,7^> =2(^,7^}, (6.12) 
1(^,7^-^,7^)1 <2||V'|||M'II 

<c||^||||9 r >|| + c||v-||||(i + ^)-'/Vll- 

From lemma l4~5l the energy controls the weighted term arising from 7, 

\\(l + rl)~ a/2 ^\\ < 11(1 + (r*) 2 )-i VII <CE[^]. (6.13) 

This proves that 7 is energy bounded. Taking a = 2 > 1 for simplicity, integrating the Heisenberg 
relation, and using (|6.1[) and (|6 . 13[) . we have that for a solution to the linear wave equation, 

(<f>,[H,'y]<f>)dT = 2(<f>M)\ t o (6-14) 
Wil + rlr^fdr^CE^,^, 
The nonlinear contribution from the Heisenberg formula is 

I - (4>nfF'<t>) + (/f'0,7$| < E[<t>,^\\fF'<t>\\ 

Using this result and estimate (|6.1|) . we can conclude 



CEfai] > 1 1 1|(1 + r 2 )-V|| 2 + ||(1 + rtytdrM* - E[4>,m\fF'tt>\\dT. (6.15) 
Jo 

This proves equation (|1.6j) in the linear case. 
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We now rerun the argument with 

r6(r.-(aoo).) 



: dr, (6.16) 



(i + H> 

In the term — <?V/ = ~gV + gl 2 V[, we treat —gV as an error term to be controlled and —gl 2 V[ as a good 
term. The term— gV' is negative between (ao)* and (aoo)*, which is a compact region, so —gV'ip 2 can be 
controlled by (1 + r 2 )^ 1 ^ 2 in the previous estimate. We can now use —gl 2 V' L as in the previous estimate to 
get positivity and to get an extra —gl 2 Vf term, so that we have (|6.15[) with an extra term. 



CE[M]> I (11(1 +rl)- 1 4>\\ 2 + ||(1 + rl)-id r M 2 

+ \\{- 9 V L ) 1 iV w cp\\ 2 - E^^WfF'^dr. 



7 Angular refinements 

We use two methods to gain additional angular regularity: angular modulation and phase space localization. 
Angular modulation involves rescaling the radial variable using the angular Laplacian operator. For phase 
space localization, we localize the observables in regions where the radial derivative is comparable to powers 
of the angular derivative. The angular modulation alone allows us to regain control over | derivatives in 
L 2 . From the local decay result, we already control a full derivative in the radial direction. The phase space 
localization allows is to work in "bands" where powers of the angular derivative are related to the radial 
derivative. This allows us to translate the radial derivative control in to control of an angular derivative 
with the loss of 6 derivatives from the width of the "bands" . There is also a e loss from the shape of the 
localization function. We combine these as a e loss of angular derivatives. 

To compute commutators involving the phase space localization, we use a commutator expansion lemma. 
To use this we need to introduce a notion of lower order, which we will take to mean local terms involving 
fewer powers of the angular derivative. 

Definition 7.1. The operator L is defined by the spectral theorem as 

L = Vl-Asnj (7.1) 

The phase space variables are defined by 

y —T m r 

Definition 7.2. The phrase "lower order terms" will refer to terms of the form 

\\L? +n X uf, 

and terms which are controlled by the local decay estimate such as 

llx^ll 2 , WxdrM\ HVauxV'll 2 , 

where \ is a compactly supported function, and x is a compactly supported function which vanishes at least 
linearly at r* = 0. These are denoted 

ogil^WII 2 ^ 1 ) 

The estimate for the phase space observable involves localization in the phase space variables. For this 
purpose, the following notation is used when the localizing function is a characteristic function. 
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Definition 7.3. For 5 > 0, the sharp, near localization and sharp, interval localization for £ n are 

%i<i =x([o,i],0 
*i-'<i«i<i =x([r s A],0- 

From these, 5>|£ |<i and ^i-a<|^ |<i can &e defined by the spectral theorem. 

We now define the phase space observable, T, and its pieces, T n , m - To do this, we will first need to 
introduce a more regular phase space localization function. Recall that Xa is defined to be a positive, 
smooth, compactly supported function which dominates the trapping terms. Essentially, it localizes near 
r = a. 

Definition 7.4. For fixed values of S > and e > 0, the following definitions are made. For simplicity, we 
will take 6 to be the inverse of a positive, even integer. 

The angularly modulated Morawetz multiplier is defined in terms of g from (|6.16[) by 

g L m =g(L m r lf ) = g(x m ), 

7l- ^^(9L"-d r , +d r ,g L m). 

The first smooth, near phase space localization functions is defined by 

$ a , e {x) ={l + x 2 )-^. 

The partial phase space ohservables are defined for < n < m < ^ by 

The phase space observable is defined for 5 — e by 

26 z 

j=o 

With the local decay observable, the phase space observable majorates L 1 " £ . Alone the commutator of 
the phase space observable dominates 1 — e angular derivatives localized for | — id Tr \ < L, up to lower order 
terms. The remaining regions of phase space and the lower order terms can be controlled by the local decay 
observable. For our later arguments, the estimate on the phase space observable alone will be necessary. 

Lemma 7.5. For ip <G §, S > 0, there is a Cr > such that 

(V, [H,T + Cn#> >C{\\L^ x n 2 + 11(1 + a-VH 2 ) (7.2) 

Proof. Each partial phase space observable majorates a power of L in a region of phase space. The partial 
phase space observables appearing in the phase space observable are chosen so that, together, they cover the 
entire phase space, up to lower order terms. Since the lower order terms are controlled by the local decay 
observable, it is sufficient to add a sufficiently large multiple of the local decay observable to get a strictly 
positive commutator. 

The estimates on the partial phase space observables are reproduced from [5] in appendix [SJ Lemmas 
IA.1I and IA.2I sav that the partial phase space observables with m = n dominate L 1_<5_ f in strips of width 
L 5 in £ n , and that the partial phase space observables with n= \ dominate ~ for |£ n | < 1. 

ftMtf,r ni ity) >c\\L^^^ n ^ lXa n 2 -o^L-^x^f.L 1 ) 



19 



Summing over the partial phase space observables, r n>TO , in the phase space observable, T, gives an 
estimate in many regions of phase space. The highest powers of L appearing in the lower order terms are 
for n = \ - 26. 

(i>,[H,r]i>) ycdiL 1 ^^, 2sl<lXa iP\\ 2 + E H^-'^-^i^Krx^m + odli 1 - 2 ^!!^ 1 ) 

5 0=0 

The regions |£i_ 2(5 | < 1 and L _<5 < |^| < 1 correspond to | - id r J < L? +2S = L 1_ ^~ 2<5 ' and 
£i-U+i)<5 < | — i,9 r> | < L x ~i s . Since j varies from to — 2, the regions covered, taken together, cover all 
of the region | — id rr \ < L. 

(iM#,iW >c\\L 1 ~^^_ tL - ldrA < lXa n 2 + o{\\L 1 - 2S xi>lL 1 ) 

The region | — id rr \ > L is covered by the local decay observable. 

(V, [F,7]V>> >IIXa^^|| 2 + ||x^H 2 

>ll^,Xa^H 2 

>||L(1 - $|- tL -ia r j<i)XaV'l| 2 - 
Putting the estimates in the different regions of phase space together, we have that 

(V, [h, r + ci\i>) ^cWL'-^x^f + o{\\i}- 25 x n,L l ). 

Since 1 — 4£ > 1 — 25, the first term dominates the lower order terms for sufficiently large L. Since the 
local decay observable majorates ||(1 + r 2 )" 1 ^ 2 > CHxV'l! 2 with no remainder, by interpolation, there's a 
constant Cr for which 

[H, T + Crim ^CiWL^x^f + ||(1 + r 2 )- VII 2 )- 

□ 

This result can now be integrated, using the Heiscnbcrg relation to conclude the L 2 space-time intcgra- 
bility result for 1 — e angular derivatives. 

Proposition 7.6. If cp £ § is a solution to the linear wave equation (|2.1[) and e > 

/>oo 



Proof. As in the proof of equation (|6.14[) . we apply the Heisenberg relation and integrate the commutator 
estimate 2[) . To do this, we need to show T + Cr7 is energy bounded. From (|6.12[) , we already have a 
bound on the term coming from 7 at the initial and final time. Since T is a finite sum of T ntTn terms, it is 
sufficient to estimate each r njTn term individually. From lemma [A.3[ for any function tp G §, 

||r„, m ^|| 2 < c n , m E[^] 

We can now integrate the Heisenberg relation with A = T + Crj- 

j t ((l A4>) - (0, A4>)) =(4>, [H, A]<j>) (7.3) 
2(4, (r + C T1 )4>) \\ > f WL^x^W 2 + ||(1 + rl)- V|| 2 dr. 



Estimating the right hand side by the energy, using the Hardy type estimate (|4.5j) . and taking the limit as 
t — > 00. 

/>oo 

CE[</>,4>]> / ll^-^X^lP + IKl + r 2 )- 1 ^! 2 ^. 
Jo 

We take e = (3/2)6. □ 
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8 Temporal Refinements 



We now apply the method from [6j to gain additional control in time. Instead of introducing an energy 
bounded observable, we introduce an observable which is controlled by the conformal energy, and use a boot 
strap argument to control the growth of the conformal energy. 

We begin by introducing a localization deep inside the light cone and another concept of lower order 
based on this localization. 

Definition 8.1. Let X[o,i] '■ ^ — > [0,1] be a smooth, even function which is identically one in [—1,1], 
identically zero outside [—2,2] and decreases smoothly in [1,2]. If no argument is given, it is assumed that 

X = X[o,i](^). (8.1) 

The constant 10 in this definition is not particularly relevant. With this, or other choices of constant greater 
than one, x wiJZ be called the light cone localization. 

The phrase "lower order in the light cone" will refer to terms of the form, or bounded by, 

E[x4>,xi>], 

and will be denoted 

0(X). 

Note that, because wc allow different constants to appear in the definition of \, and 
„ ,„ ,10r* 10 w ,10r C „ 5r 

we write d Tt x = (1 + As long as only finitely many derivatives of x are taken, we can control x an d 

its derivatives by x taken with a different constant in the argument. Each time a derivative is taken, because 
the support of X[o,i] is in \x\ < 2, the constant in the argument decrease by a factor of 2. By starting with 
a sufficiently large constant in equation (|8.ip . we can ensure that, even after dominating derivatives, the 
constant in \ wm be larger than 2, and the support of x wu l be strictly inside the light cone. In particular, 
we will never take more than two derivatives, so the choice of constant 10 in (|8.1|) is sufficient. Were we to 
require control on more derivatives of x, we would simply take a larger constant in (18. ip . 

By the Poincare type estimate in lemma 14.51 we can also control weighted norms of <f> with light cone 
localization. 



2 



1 



cPfi <CE[ip,ip] 



(xVO d> <CE[x^,m = o(x). 



1 + ri 

This also means that the order of x localization and derivatives does not matter. 

<MxV>) =X[o,i](:p^R,V' + 7^X[o,i](yq7pV> 

<x [0 ,i ] { Tr Jd r ^ + — X[0 , 1] {—^ 

I 10r * \a i , C - t hr * \ i 
<X [ o, 1 ]( r ^)9 r ^ + ^yrX[o,i]( TTi )^ 

\\d r ,(m\ 2 =0(X)- 

Using this light cone localization, we introduce the temporal phase space observable. The first part is 
(1 +i) times the phase space observable, and the second is (1 + 1) times the local decay observable localized 
inside the light cone. There is no need to localize the phase space observable inside the light cone since it 
already contains the localization Xa, which is localized near the photon sphere. 
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Definition 8.2. The temporal phase space observable is 



f =(l + t)r + C r (l+t)x(^)7X(; 10: 



•1 + t' !/ ^l + t' 

Note that T is anti self adjoint since 1 + t and x are bounded functions at each time, and T and 7 are 
anti-self-adjoint. 

The temporal phase space observable majorates 1 — e angular derivatives and ^ factors of t; that is, the 
commutator dominates 2 — 2e angular derivatives and one factors of t in expectation value. To show this, 
we must also include the contribution from the time derivative. 

Lemma 8.3. For ip 6 §, and e > 0, 

(V, [h, r» + (V', U) - (i>, >(i + t)c(||L 1 - £ x ^ll 2 + ll(i + ^)~ 1 x^ll 2 ) + o(x) 

Proof. This lemma says that f majorates t'L 1-6 . Since f is essentially tr, this should be expected. In 
essence, the left hand side is t times the commutator with V, so that the result follows from proposition [731 
There are corrections from the commutator with the light cone localization in T and from the time derivative 
of f . These are shown to be lower order inside the light cone. 

The commutator can be broken into four pieces using the linearity of the commutator and the Leibniz 
rule for commutators. 

[H, f ] =(1 + t)[H, T] + C r (l + t)x[H, 7 ]x + C r (l+ t)[H, xM + W + t)xi[H, x] (8.2) 

We now show that the first two terms provide the bound that we want, that the three remaining terms are 
lower order in the light cone, and that the terms involving T are also lower order in the light cone. Combining 
these three steps proves the desired result. 

Step 1: The dominant terms in the commutator The contribution from the first two terms 
can be found using proposition 17.51 



(V, (1 + t)[H, T] + C r (l + t)x[H, 7 ]xV'> >(!+*) (qiL 1 -^,^-^ |< lX o«|| 2 + OQIL 1 ' 2 ^^ L 1 ) 

+ C T (^,x[H, 1 ]xip}) 
Since (1 + r^)^ 1 ^ > Cxa, the interpolation in proposition 17.51 still holds, and 



/ 1 0r 

Oj, (1 + t)[H, T] +C r (l + t)x[H,l\xm >(1 + t)C UL^x^t + 11(1 + r 2 X l xi^f-^A 

^(l+WWL^x^f- 

Step 2: Remaining terms in the commutator We now consider the last two terms in equation 
The fourth term is, up to signs, the adjoint of the third. When the derivative in the commutator 
lands on Xi a negative power of (1 + t) is introduced, so that the commutator is lower order in time. It is 
also localized deep in the light cone. 

(1 + t)[H,x] -(1 + t) (-2d r MM^- t )) + (dlx(^- t )) 

, n1 .,A0n, 100 



1 + t' 1 + t 
105 w 

ix =xi + YTl x 

Since 7% and the adjoint of {1+t) [H, x] both look like the sum of terms involving one or zero radial derivatives 
times bounded functions and light cone localized functions, the remaining commutator terms are lower order 
inside the light cone. 

(v, (1 + mH, xhm + (tp, (1 + t)xi[H, m =o{x). 
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Step 3: Controlling the T terms Since T is anti-self-adjoint, it is sufficient to consider 
We first compute f . 

- lOr 

r = r + c r xix - C r (i + t) T rrhW'rt + xix')- 

{i + ty 

We will deal with each piece separately. 

Since Xa is localized inside a fixed compact region, and x is identically one in a larger region for alH > 1, 
X localization can be added to T for free. 

r = xr 

Therefore, 

<llx^lll|iYH 

From lemma |AT3| 1 1 r-*/> 1 1 2 is controlled by the energy of Xa">P, so 

(hry) <\\m\E\x a ^ 

<\\m 2 +Eim 

<o(x) 

Clearly, the XIX terms are lower order inside the light cone. 

(ip,xixi>) <\\x^\\\hx^\\ =o(x) 

The remaining terms, from when the time derivative falls on x, are dealt with similarly. The derivative 
gives an extra factor of + t). This is bounded inside the light cone. 

lf)r Ir I Ir I 

□ 

In our proof, we bound terms both by the energy and the conformal charge. Energy boundedness was 
central for proving equation (|6.14|) . and we require an analogue combining E and Ec- For ip £ S, Q(t) 
depending on t and ip is conformally bounded if 

Q(t) <min{(l + i)£;[^,^,(l + i)- 1 i? c [^,V;]}. 

We say an operator A is conformally bounded if ||^4'(/;|| 2 is conformally bounded. 

We now turn to proving the main estimate, that the conformal charge of <f> and the energy of L £ <j) remain 
bounded for the linear wave equation or for the nonlinear wave equation with small initial data. In the linear 
case, we can apply L e to the wave equation to show that the energy of L £ (f> is conserved, and then apply 
the angular refinement argument to control Lcf> in L 2 . This controls the growth of the conformal charge. In 
the nonlinear case, we can still use the energy of L £ <f) to control the growth of the conformal charge. The 
conformal charge can then be used to control weighted L q norms, which in turn control the growth of the 
energy of L e (j). Through a bootstrap, for small data, the energy of L e <p and the conformal charge control 
each other. 

Proposition 8.4. If e > and <j> £§ is a solution to the linear wave equation (|2.1|) . 

EcW),kt)\ <cE C [m,m]+cE[L £ m,L £ m] 

For 8/3 < p < 3 and sufficiently small initial data, in the norms appearing on the right hand side, the 
same estimate is true if (f> £ S is a solution to the non-linear wave equation (|2.2I) . 
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Proof. If it weren't for the loss of e angular derivatives on the commutator estimate of lemma 18.31 we'd 
expect to be able to close the estimate on the growth of the conformal charge using conformal boundedness 
of the operator T. Instead, we will have to resort to a more complicated notion of L e -conformal boundedness. 
Our method follows from the following six estimates, in terms of a constant TV to be chosen later, 

E[L 18e (f>, L 18e <p] (t) <E[L 18e (f>, L 18e 0] (0) Persistence of regularity (8.3) 

+ f E[L 18 ^, L ls ^]i\\L^f(n)F'(\cb\ 2 m\dT 
Jo 

E c [(f>,(j)](t) <E c [<t>,<j>](0) +C / t\\LxcM\ 2 oIt Conformal estimate (8.4) 

Jo 

+ C [ r [ (2f + r,f')F(\^\ 2 )d 3 f idr 
Jo Jm 

[ T\\L Xa <p\\ 2 dT <CN 2 sup £;[i 18e </-,L 18e 0](r) Temporal refinement (8.5) 

JO 0<T<t 

+ CN~i sup E c (t)+C f |(/F'0,L 2e f0)|dr 

0<r<t Jo 

[ E[L 18e <f>,L 18e <f>}?\\L 18£ f(r^)F'(\(t)\ 2 )(t)\\dT NL persistence error (8.6) 

Jo 



<C sup E[L 1S£ cl),L 1S£ (j)}(T) sup E c [4>,$\{t r 

0<T<t 0<T<t 



[ t [ (2/ + nf')F(\<f>\ 2 )d 3 fidT NL conformal error (8.7) 

Jo Jm 



<CN^E + CN^r 1 SU p E c {t)^ 

0<T<t 



[ (fF'(f),L 2s T<f>)dr <C sup E[L 18e <j>,L 18e <j)](T)? sup E c [(/>, <j>] (r) § NL Morawetz error (8.8) 

JO 0<T<t 0<T<t 



The conformal estimate (|8.4[) is already known from 

In the linear case, the persistence of regularity is trivial since L 18s commutes with the d'Alembertian, so 
the E[L 18e (f>, L <j>] is bounded by its initial value. The non-linear estimates l|8.6p - (|8.8[) are not needed in the 
linear case. We note the loss of angular regularity comes entirely from (|8.5[) as a consequence of the probable 
lack of angular derivative decay near the photon sphere. There is a lose of 2e in the commutator estimate, 
and, in distributing, this between the energy and the conformal charge, we put 18 factors on E[L 18e (f>, L 18e cj)] 
so that we can get none on Ec and get a factor of ./V -1 / 4 . 

Step 1: Persistence of regularity, (|8.3[) - For this step, a will be an index of angular regularity, 
and not the value of r corresponding to the peak of the effective potential or to the photon sphere. We apply 
the operator L a to the wave equation, 

d 2 (L a ^) = dl(L a cf>) - y(L a 0) - V L (-A w ){L a cf>) - L a (/F'0), 

and consider the densities and relations with general nonlincarity J\f — L a (fF'(f>), 

e [v, w] =7j( w2 + (d r ,v) 2 + Vv 2 + V L \W w v\ 2 ), Po*[ v , w ] =wd rt v, p QuJ =wV w v. 

If if) solves the general nonlinear wave equation, 

2 (iP) = d 2 t (iP) - V(if)) - V L (-A S2 )(i))-Af, 
then there is a relation for eo, and a formula for the time derivative of its spatial integral. 

d d f f 

d t eo[tp,ip] =d r ,po*[ip,^] + Van ■ pou[if>,ip] - —E [if),ip](i) =— e [if>, ip}d 3 fi = - / ipj\f. 

at at J m J OT 
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Taking ip = L a (f), 

E[L a d Pl L a ^](t)=E^ 1 i,]{t)= f f L a <pL a fF'(\^d 3 ^dT + E[L a ^L a ^}(Q) 

Jo Jm 

< [ \\L a 4\\L a fF\\<j>\ 2 )ndT + E[L a ct>,L a ti(0). 
Jo 

We now take a = 18e. 

Step 2: Temporal refinement, (|8.5p . We consider the operator L 2e f and apply the Heisenberg 
identity 

IdlL^fcf,) - (cb,L 2 ^H)) ={4>, \H^t]4>) + (0,L 2e f0) - (cf>,L^H) 

+ (fF'4>, L 2e T(j)) - (0, L^tfF'cp) (8.9) 

Ignoring factors of L e , our method is to show that the term on the left, under the time derivative is 
conformally bounded and that the linear term on the right dominate tL 2 x a in expectation value, so that the 
time integral of the expectation value of tL 2 is conformally bounded. We now show this with the additional 
factors of L e as required. The last two terms in equation (|8.9[) are not present in the linear case. 

We begin by bounding the skew expectation at the end points. Since f = (1 + t)(T + CrXlXii an d 
f = fx, we can estimate the norm of f, ||r>|| 2 = ||fx</>|| 2 < (1 + t)E[x4>,x4>] < E c [^,4>]. Since L 2e f is 
anti-self-adjoint and additional light cone localization can be applied to T, 

|(0,L 2e f» - (0,L 2£ f»| =2|(L 2e x0,i>)| 

<C£[L 2e 0,L 2e #£ c [</>,#■ (8.10) 

We proceed temporarily using the notation 0(L £ x) to denote 0(x) with respect to L e (j> instead of </>. 
From lemma [B~3l applied to L £ (f>, since L e commutes with T and H, 

(0, [if, L 2e f ]4>) + (0, L 26 t <f>) - (0, L 2e t 0) >C(4>, tL 2 Xa & + 0(L e x) (8.11) 
The nonlinear terms are 

|(/FV,£ 2£ f0) - (<f>,L 2 *rfF'<f>)\ <2|(f0,xL 2e /F'(|0| 2 )0)|. (8.12) 
We integrate the Heisenberg relation (|8.9[) and substitute (|8.10p - (|8.12p . 

CElZFhZFfiiEcfaifi > f \^TL 2 Xa <f>)+0(L £ x)dT 

Jo 

-C\( (fF'^L 2s f^)\dT 
Jo 

Solving for the time integral of tL 2 , we get 

/ (&TL 2 Xa 0)dT <C sup £[L 2 ^,L 2e 0](r)s sup E c {4>,^]{r)^ 

JO 0<T<t 0<T<t 

+ f 0(L'x)d T 
Jo 

+ C\ f\fF'^L 2s U)dT\. (8.13) 
Jo 

To prove (|8.5p . it remains to bound the integral on the right. We decompose the integral onto harmonics 
and then decompose into small and large time relative to the harmonic parameter. This allows us to make 



25 



estimates based on the supremum of norms of <j> rather than the sum of suprema of the projections of </>. We 
recall that P; denotes the projection onto the I th harmonic. We use A 2 to refer to the eigenvalue of 1 — Agrr- 

0(L e x)dT <C^A 2e / E[P^,P^]dT 
Jo l Jo 

^ Jnx 8 ? 

<cViV 2 A 18E / £;[P i 0,P i( /)](T)(l+T)- 2 dr + cViV-W EclP^P^jT-idT 
x Jo x Jn\ 8 ? 

<CN 2 f E[L 1Ss <j>,L 18e <j>}(T)(l +T)- 2 dr + CN-i f Ec\^\r~^dT 
Jo Jo 

<CN 2 sup E[L 18e 4>,L 18e <P](t) [ (l + T)- 2 dr + N-i sup £ c [0,0](r) [ (l+r)-*dr 

0<T<t Jo 0<T<t Jo 

<CN 2 sup E[L 18s 4>,L 1Ss 4>](t) +CN~* sup E c [4>,<j>](T) (8.14) 

0<T<t 0<T<t 

The contribution for t < NX 8e can be reduced to NE[L 10£ (f>, L 10e 4>] in the linear case, but there's no 
major gain in the argument. From (]8 . 13[) and (|8.14[) , 

ft 



f {<P,L 2 Xa 4>)dT <C sup E{L 2 ^,L 2s <j>}(T)? sup E c [^4>]{t)^ 

JO 0<T<t 0<T<t 



+ CN 2 sup E[L 18e ^,L 1S£ <j>](T) + CN~* sup E c [^,4>](t) 

0<T<t 0<T<t 
t 

2ei 



+ C\ / (/F'0,L 2e f</»)dr| 
Jo 

ft 

<CN 2 sup E[L 18e 4>, L 18e (p](T) + CN~^ sup E c [(f>, 4>]{t) + C\ / (fF'<j), L^T^dr] 

(8.15) 

Step 3: NL persistence of regularity, (|8.6[) . We continue working with an arbitrary angular 
regularity exponent a. We begin by stating a fractional chain rule estimate in R™ Proposition 5.1, pll2]. 
In E", with s 6 (0, 1), l/q = l/qi + l/q 2 , and x i-> [x^x a C 1 function, 

Illwr^lU-™ <C||Mf- 1 || £ « 1 ||u|| H ..« a . (8.16) 

Here H s ' q denotes the space of functions with s derivatives in L q , where s derivatives are defined by the 
Fourier multiplier (1 + |£| 2 ) s / 2 . By taking a partition of unity, we have the same result on the sphere. Since 
L corresponds to the Bessel potential (1 + l^) 2 ) 1 / 2 = (1 — A) 1 / 2 , taking qo — 2, q\ — 3, and qi = 6, and then 
integrating in space against the weight / wc have 

||L 18E /(r*)F'(M 2 M| ^H/F-^r^UII^^He. 

Under condition [6l the L 6 norm of L 18e </> can be estimated by the Sobolev estimate (|4.1|) with the non- 
derivative terms estimated by the Hardy inequality (|4. T|) . For 2 < 3(p— 1) < 6, the second can be estimated 
by the the conformal-Sobolev estimates (|4. 3[) . 



|L 18e /(r*)F'(|0| 2 )0|| KElL^^L^^Ec^JrHl + t)-^- 
<E[L 18e 4>,L W£ 4>]E c [4>, fif-^l+tyP+i 
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This can now be integrated in time if p > | , 

' ||L 18£ /(r.)F'(0 2 )0||dr < f E[L 18 ^, L 18 ^]£ c [0, 4>}^{1 + t)~*+Ut 

Jo 

< sup E{L 1S£ <j),L 1S£ <}>}(T) sup EelMW* 1 

0<T<t 0<T<* 

It is at this stage that we have used the hypothesis 8/3 < p. 

Step 4: NL Conformal, (|8.7|) . Since F(|</>| 2 ) is positive, and 2/ + r*/' is bounded by a compactly 
supported function, and hence by / itself, 

(2/ + r*f')F(\<l>\ 2 )d 3 n < C [ /F(|^| 2 )dV < CE. 

For this argument, we use x to denote compactly supported, positive functions, which may vary from 
line to line, with C or B for constants or bounded operators. Let \ be a function which dominates 2/ + ?%/', 



(2f + r*f)FM 2 )d 3 v< / xWIVm 
ot Jm 

Since F is a power semi-linearity, up to constants -F(|<?j| 2 ) = CF'(|0| 2 )|</>| 2 , so 

X F(H 2 )dV>c||x0||||xnH 2 )0ll 



/ an 

These can be estimated by (|4.4[) and (|4.3p . 

Note that more powers of t can be controlled than usual because x</> is compactly supported, so that Hx^ll 
and \\d rt (x4>)\\ both decay like t~ 1 Ec[4>,<ft], instead of being merely bounded like and ||o*r.0||- 
Integrating in time, we have 

t t 

[ t[ (2f + rJ')F(\<j)\ 2 )d 3 [idT <C [ Erdr + [ i E^T- p dr 
Jo Jm Ji Jn? 

KCN^E + CN^r 1 sup Ec{t)^ 

0<T<t 

Note that, if we had merely required that the nonlinear trapping terms were merely bounded by /, then 
we would only control \\<fi\\ and ||<9 r „</>|| by Ec[4>, </>]. In this case, we would lose two factors of r in the estimate, 
and have t 2 ~ p in the integrand. Since the Sobolev estimate requires p < 3, this would not be integrable. 

Step 5: NL Morawetz, ()8.8|) . The nonlinear terms can be estimated by using the Sobolev estimate 
(|4.1[) and conformal-Sobolev estimate (|4.3[) as in the persistence of regularity argument. 

<ciif0iniL 2£ /i0rvii 

We now estimate ||f 4>\\ by 

\\Y<j>\\<CE c m,M<CEc[&4>\, 

and ||L 2e /|0| p_1 (/)|| the same way as in the non linear persistence of regularity argument, and integrate all 
of this in time. 

\\L a f{r*)F\4> 2 )<t>\\ <E[L*tj>,L*$Ec\M}^{l+t)^, 
r|(/F'(|0| 2 )0,L 2 -f0)-(0,L 2s f/F'(|0| 2 )0)|dr < sup i?[L 18 ^, L 18 ^](r)^ sup Ec[M](t)S 

Jo 0<T<t 0<r<t 
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Step 6: Closing the boot strap. We summarize the results of estimates (|8.3[) - (|8.7[) as 

E[L 1S£ 0,L 1SE <f>}(t) <E[L ls ^,L 18s 4>}{0) + C sup E[L 18e <f>, L 18e 0] (r) sup 0](r)^, 

0<T<i 0<T<t 

Ec[MW) <Ec[(f>,j>](0) 

+ CN 2 sup E[L 18£ (I>,L 18£ ^}(t) + CN-t: sup £ c [<^](t) 

0<r<t 0<r<t 

+ C sup S[L 18e 0,L 18£ 0](r)5 sup E c [<t>,<j>]$ 

0<T<t 0<T<t 

+ CiV3S[^,^] + sup Ecfiir) 2 ^ 

0<T<t 

In terms of the norm 

((f>(t),(t>(t)) 2 =N?E+ sup E c [<f>,<j>]{T)+N* sup £[L 18E 0,L 18e 0](r), 



(8.17) 



0<T<t 



0<r<t 



these two estimates can be combined as 





< 




2 i 

+ CN~* 




2 i 

+ CN* 





p+l 



We now take N sufficiently large, so that CN 1 / i < i. 



<c 



(0(0), 0(0)) 



p+i 



For sufficiently small initial data in the 



norm, which means in both the 



18ej, r.l8e. 



and £"c [ 



norms, the 



are also uniformly bounded. 



norm is uniformly bounded. As a consequence, the 2£[X 



- 18e 



and Be [ 



norms 
□ 



From the boundedness of the conformal charge, the Sobolev estimate (|4.1| and the control on the L 2 
norm by conformal charge in equation (|4.4|) in the proof of (|4.2p , the L 6 norm decays 



F L 2 |0| 6 dV<£g(i + O 



an 



2\-2 



Because of the time reversibility of the wave equation, the same estimate holds as t — > — oo. 
For the wave equation on the warp product manifolds considered, this means that 



an 



an 



<( l^NW v£\<t>\»d A n) 



an 



an 



<E C (l + i : 



2\-l 



Since t~ 2 is time integrable and Eq is bounded, the spatial integral of |</>| 4 is time intcgrable. Thus, |</>| 4 is 
integrablc in space-time. 

For the wave equation on the Schwarzschild solution, the situation is a little more complicated. We 
consider the space-time L 4 norm directly, and recall that the natural measure has an additional factor of 
(l-2M/r). 



\(j>\ A r 2 drd 2 wdt 



ooJm 



(1 - 2M/r)r~ 2 \(j)\ i d i ^dt 



oo JWl 



< 



< 



an 



E C m),Ht)?o-+t 2 )~ 1 dt 
< c(E C [m,<p(o)]+E[L E m,L E m}) 2 - 
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8.1 Summary of estimates 

We now summarize the key estimates for linear solutions. These hold both for wave equation on the 
Schwarzschild solution and on Ricmannian manifolds under the conditions in section [2] 

J \\{l + rl)- l ^\\ 2 dT<CE (8.18) 

I V L \Vw<P\ 2 d 3 ndT ^CEci^il+t 2 )- 1 (8.19) 
Jm 

I x(n)(|0| 2 + |0'| 2 + |0| 2 + |V 2IJ ^| 2 )d 3 A1 <Ci?cW(l+< 2 )- 1 (8.20) 

f VM 6 d 3 t-L<CE e (t) 3 (l + t 2 )- 2 (8-21) 

where x( r *) refers to any compactly supported function and 

E=E[<t>(0),j>(0)] (8.22) 
E e (t) <C(Ec[<f>(0),j>{0)] +E[L^(Q),L e m\) (8.23) 

For small-data, non- linear problems, there is a correction term in the first estimate and the others remain 
unchanged. 

From these, solutions to the original, geometrically defined wave equation \3<j> = satisfy the following 
space-time integral estimate, in terms of the geometrically defined 4-volume. 

J |0| 4 d 4 vol<C(£ c [0(O),0(O)] +E[L e (jy(0),L e <j)(0)}) 2 (8.24) 

The region of integration is the exterior region, region I, on the Schwarzschild solution or the entire space-time 
K x 9Jt on a Riemannian manifold. On a Riemmanian manifold, for the semi-linear equation \3<p = \(f)\ p ~' 1 (f) 
with small initial data and p G (8/3, 3), the same estimate holds. 



A Estimates on the partial phase space observables 

We present three results from our previous work [5] . These are phase space localized estimate which we use 
to build the estimates involving the phase space observable T and reduce the loss of angular regularity at 
the photon sphere to e derivatives. The first two estimates state that the partial phase space observables 
majorate powers of L, and the third, that the partial phase space observables are energy bounded. For these 
recall that r il rn is defined in terms of S and e. 

Lemma A.l. For ip E §>, 5, e > 0, and n € [0, there is a constant C such that 

(1>, [H, r„ 4 ]v) >qjL lz ^$i a < lXc ^|| 2 - oiWL^xaiPlWL 1 ) 

Lemma A. 2. For G S, S, e > 0, and n G [0, |], 

{if,, [H,T n , n ]i>) >C||L 1 -4-^ i -,<| € j< lX aV'l| 2 -0(l|i i ^ i XoV'l| 2 ,i 1 ) 
Lemma A. 3. For ^GS and n < m < 1/2, T n ^ m is energy bounded: 

Cnm || r \i,m ^|| 2 <£[V,# 

In section [3 we define the operators L, x m , and £ n ; the concept of lower order, 0(\\L~ ^~ Xa^\? i L 1 ); 
the phase space localizations, $ a ,e(^ n ), |<i, v E , i-5<|^ n |<i; and the partial phase space observables, T n . m . 

Lemmas I A. Ill A.3l are essentially phase space localized versions of the estimates on 7 from section [6l Since 
the partial phase space observables are phase space localized versions of 7, with extra powers of L, the results 
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would follow immediately if phase space localization could be commuted through spatial localization freely. 
The essence of the proofs in this section is that commuting the phase space localization terms only generates 
lower order terms. We state the general commutator expansion theorem of Sigal-Soffer |29| and specialize 
to the case [Fi(x m ), F 2 (£ n )} = L m+n ~ 1 B. 

We use the notation to denote the Fourier transform. The kth commutator, Ad A (B), is the quadratic 
form (Ad^ _1 (-B)^, Aip) — (Aip, Ad^ _1 (B)ip) . If Ad A extends to a bounded operator, the same notation refers 
to this operator. 

Theorem A. 4 (Commutator Expansion Theorem). If n > is an integer, A is a self-adjoint operator, F\ 
is a self-adjoint operator for which for 1 < k < n, Ad A (Fi) extends to a bounded operator, and F2{x) is a 
smooth function satisfying W^FlF^jWi < oo, where J-[-} denotes the Fourier transform, then if [F\, F2(A)] is 
defined as a form on the domain of A n , 

n—X -. 

[F lt F 2 (A)] = ^Fl k] (A)Ad k A (F 1 ) + R n (A.l) 
fe=i 

in the form sense with the remainder R n satisfying 

\\Rn\\ <qi^i n] ]||iP<Cl(*i)ll 

Consequently, [F\, F2(A)] defines an operator on the domain of A n ~ x . 

We now specialize to the case that F\ is a function of x m and A = £ n . The commutator [x m , £ n ] can be 
computed explicitly as iL m+n_1 which has a negative power of L for n < m < i. 

Lemma A. 5. Let k be a positive integer, F\(x) be a smooth function with, for 1 < j < k, ll-Fi^lloo < °°; 
and F2(x) be a smooth function with H-^l-F^jlli < oo, then there is a constant Ck depending only on k, such 
that 

fe-i 

Lfc (l-™-n)[ Fl(xm ) 5 F 2{U] = £ p L (k-m-m-n) F [0] ^pti] ( £j + ^ 

3=1 

\\R k \\ ^Cfcii^Hooii^ 1 ]!!!. 

Proof. This follows from making a harmonic decomposition, where L = A, and noting the adjoint is 

Ad{ n (F, (x m )) =iA" _1 9 r .Ad^~ 1 (Fi (x m )) 
=p \3( n -V (d r ,Y F^r*) 
=i j \i( m+n - 1 ') (A m r* ) , 

and applying the commutator expansion theorem. □ 

The operators gL m , &a,e{€ n ) an< ^ |<i an< ^ defined in section [71 We also introduce the following 
functions. 



{l + b\x\) 21 1 + x 2 



X l( x ) = 7i , UZi\2 ' ^l(- T ) = i , „2 ' x t(- t ) = x 



The operators Aj(x m ), A|(x m ), X-f(x m ), <5>£, je (£ n ), an d $c,e(£ n ) are defined by the spectral theorem. Two 
important properties of these functions are 

X^xf+X^xf >C, (A.2) 
*«. e (0 >C$ a , e (0 (A.3) 
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In applying the order reduction lemma, it is useful to note that, if Fi is a C 1 function with bounded 
derivative, and if $(£) is one of xi, $ 0ie (£), 4>& iC (£), or <j> Cje (£), then 

[*!(*), *(0] =L n+m ~ 1 B 1 

Bi <q|P?||oo, (A.4) 

where the notation B is used to denote an arbitrary bounded operator which commutes with L, the same 
way that C is use to denote an arbitrary constant. In addition, if F\ has a bounded second derivative, then 
for*(0=£*„ ie (0. 

[^(x m ), *(€„)] =L n+m ~ 1 B 
USHCCKIU + UFflU). (A.5) 

In particular, equations (|A.4[) and (|A.5[) hold for i*\ being any Schwartz class function, Xj^x), or X|-(x). 
These estimates follow from the order reduction lemma and the fact that the Fourier transform of $' is in 
L 1 because the high degree of regularity for <E> translates into decay properties for the Fourier transform. For 
= x§ a ^{x), has Fourier transform in L 1 and $' in L°°. 

We now prove lemma IA.3I The method of commuting localization to get lower order terms will be 
typical of later arguments. We will also use the notation Qi < Q2 to denote there is a constant C such that 
Qi < CQ 2 . 

Proof of lemma YA.S\ We start by dropping a spatial localization and expanding 7^™ 

HiWII 2 =Wx a ^aAUL n ^iL^aAUx a M\ 2 

We rewrite the first term in terms of £ n and drop a bounded operator in the second, leaving the bound 

We commute phase space localizations & a A£n)£n an d 9L m m the first term and write the commutator at 
the end of the expression, as we will typically do. The bound becomes 

WL^QL^AinfinXaXn 2 + \\L m+n -^\\ 2 + WL^^AU^, QL^aAUX^f- 

We now drop the bounded operator g^m in the first term and apply the order reduction lemma in the last 
to get the upper bound 

\\L 1 - c *aAtn) 2 S n X«'l>\\ 2 + \\L m+n -^f + \\L m + n -^ a AUX^f- 

From the decay of <& ai£ (£„), we have that 

$a, £ (o 2 e< 

Applying f 1 ~ e g e < \ f\ + \g\, and interpolating L m+n between L and 1, gives 

!|r, M „VII 2 <\\Lx a M\ 2 + Wdr.XaMl 2 + \\L Xa M\ 2 + llxc^ll 2 

□ 

To prove lemmas IA.1I and IA.21 we prove preliminary bounds on the commutators involving H . These 
bounds arc localized in both x m and £ n . The key idea, as demonstrated in the previous proof, is that 
commutators arc lower order terms. 
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Lemma A.6. For ijj £ §, n < 1/2, n < m < 1/2, and e, 5 > 0, there are constants depending only on n, m, 
e, 8, such that 



(i,,[H l +H a ,r ntm ]^) >a($ a , e (^)x^,i"" e (-^,^(x m )9 r , -£ a ff£»VZ)*a,e(0*nXa^) ( A ' 7 ) 

+ 0(||L^ i Xc( VI| 2 ,i 1 )- 



Proof. We first consider the commutator with Hi. The commutator with i?2 is expanded by the Leibniz 
rule. The commutator with \a is zero. 



[Ha,r„, OT ] -Xa[^$a,e(^)]7L™i n - e $a, e (IJXa+Xa$a,e(l„)7L™i"" e [^ < I ) a,e(^)]X ( 
+ X a * a , e (l„)[^7L™]i"- e $a,e(l Il )Xa. 



The first two terms are, up to signs, adjoint. The commutators involving < &a,e(€ n ) can be expressed using 
the order reduction lemma, and the commutator with 7^™ can be expanded as in section [6] 



\^AH2nL^)\S(^aAUXa^,L n -^V^ a AU]X^}\ 

+ \($aAUx a i>, [v, lLm ]L n -^ a AUx^)\ 
<\h$aAUx a i>\\ 2 + WL^L^BxaM 2 
+ \\®aAUx^\\ 2 + \\9L^v'L n ~^ a AUx a M\ 2 

^OiWL-^x^WL 1 ). 



By an identical argument, the difference between the commutator [H 3 ,T n ^ rn \ = [(L 2 — l)Vt,r„ im ] and 
[L 2 Vi,r„ jm ] is lower order, and we will ignore this difference in the following arguments. 
The terms coming from Hi + H 3 can also be expanded using the Leibniz rule. 



[H X + H 3 ,T n , m ] =[Hi + H 3 , Xa }^aAUL n - e 7L^aAUXa 

+ X^aAtAL^'lL^aAiAWl + H 3t Xo 
+ Xa [Hl + H 3 , QaAUW^lL^aAUx, 

+ x^aAUL n ' e -/L^ [Hi + H 3 , $„, £ (£ n )]x 

+ X«$aAUL n - e [Hl + H 3 , 7L»]$a,e(€„)X. 



The first two terms are estimated first. Since H3 = [l? — 1)Vl commutes with functions of r*, the 
contribution from H 3 can be ignored. The terms are, up to signs, each other's adjoints, so it is sufficient to 
consider the second. We work in expectation value and begin by expanding 7^™ , 



<\^,X a ^aAUL n - e d r , 9L ^ a AU[-dl,X a W\ 

+ 1 (V, x a 4>aAUL n - e L m x l (x m )$ a , e (cj[-a r 2 , , X a]V> I 

<\(dr,X a i',L n -^ a AU9L^aAin)[-dl,X a W\ 

+ \(L m+n -'$ a AUx«i>,xA^)*aAU[-dl,x«m- 



(V, [H 3l r„, ro ]v) =0(111,^ x^f-L 1 ) 



(A.6) 



\(^,X a ^aAUL n - e iL^aAU[~dl,X a m 



Expanding 



[-dr^Xa] = ~d r ,x' a ~ X, 
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and applying the Cauchy-Schwartz inequality gives 

<\\d r ,xAA\ 2 + \\L n -^ a ^U9L^aAU[-dl,X a W + \\L m + n -* X ^f 
<l|£ B - e *.,e(*„)ffL»*«, 6 (^ 

BL n -^aAU9L^aAUdr t X' a i'\\ 2 

+ \\L n - f -$aAU9L^aAUX>\\ 2 + \\dr,x'JA\ 2 + IW^f + 0( \ \ L ^ X ^ II 2 , L 1 ) . 

The second, third, and fourth terms are clearly lower order. It remains to estimate the first, which we rewrite 
and then estimate by the order reduction lemma 

|<^Xc*$a,e(0£"~ e 7^a,e(O-Cx«MI<l^^ 

+ 0(\\L 1 ^ Xa ^\\ 2 ,L 1 ) 
We now apply the order reduction lemma again using <& aie (()£) to get a bound of the form 

IMX«*.,«(0£ n ~ 6 7i»*a, e (Oh^^ 

^-^AUHnx'M^o^L^x^fiL 1 ) 

a^-% n \ e x: a n 2 +o{\\L- i ^ Xa n\L i ) 

<\\LAM 2 + \\\Ux'AA\ 2 + 0{\\L-^ Xoi i>\\\L l ) 
<\\Lx'JA\ 2 + WL'^d^i'W 2 + oiWL-^x.n 2 ^ 1 ) 

The second term is clearly lower order since L™ -1 is a bounded operator. The first is also lower order since 
the localization x' a vanishes in a neighborhood of r, = 0. 
We now turn to the remaining three terms. 

[i?i + ff 3 ,r„, m ] ^xAHi + H^aAUW-'iL^aAUxo. + x^aAUL^iLAHi + H^^aAU^ 

+ X a <£aAUL n - e {Hi+H 3 , lLm ]<t> a AUx a - (A.8) 

Since Hi = —d 2 t commutes with c E > a ,e(^„), the contribution from Hi can be dropped from the commutators 
with $ a ,e(Cn)- The contribution from the Hi commutator is 

X^aAUL n - ( [Hl,lL^aAUX a = 2x a <£aAUL n - e (-dr,g' Lm d r , - ifl#»)*., e (£ n )Xa 

The contribution from the H3 commutators is significantly more complicated than the previous pieces, 
and the commutator with <fr a ,e(£ n ) must be expanded to second order, 

[H 3> * a , e (0] HKAU[Vl,U + KAU[[Vl,U,U + Rs)(-a w ). 

Since the difference between the contributions from (— A<hj) and L 2 is lower order, we will ignore it. The 
commutator with Vl(L 2 — (— Agjr)) ~~ Vl alone can be treated the same way as the commutator with Hi- 
Using the second order expansion with three sub-terms, and expanding 7^™ as the the sum of the two 
sub-terms, for each of the first two terms in equation (|A.8j) gives twelve terms, plus one more from the third 
term in (|A.8j) . for a total of thirteen. We treat the radial derivative terms in y^m as higher order, and group 
the thirteen terms in terms of formal order, with the third term in (|A.8|) treated as the same order as the 
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top part of the other two. We strip out localization in % a , and get 

[H 3 , r n>m ] =Xa(Ql + Q2 + Q3 + Qa + Qb + QsjXa 
Q X = L 2+n - f i^ e (U[VL,Udr,9L^aAU 
+ ^a,e(^„)5L".[Vt,5 r ,]$ , e (O 

0.3 = L 2 +' n -^^AU[[VL^ n iUdr,9L^aAU 
+ $a, £ 9 r , $2 )6 (£„) [[y L , £J , 

04 = i 2+n -^( - < e (0[Vb,U€»]$U*«,«(*») 

+ *a, e (Ofl'i»< 6 (€n)[[^ > €n].€J) 
Q 5 = L 2+n - e {R z d r ,g L ^ a AU + *aASn)9L>»d rm R a ) 
.1 



The top term Q\ gives the desired estimate. All the others are lower order. We must expand to second 
order so that the estimates on the remainder, which can not utilize the structure of the localization, arc 
lower order in equations (jA.llj) and (|A.12j) . 

We rearrange the terms in Q\. It is at this point that the positivity of e and the smoothness of the 
localizations in £ n are used. 

L -(2+n-e) Qi ^AUi-VlX-iL^dr.g^AU + SaAUi-VliL^g^aAU 

+ *-, e (€„)ff£»^X, e (€„)(-V£) 

We rewrite this in terms of and commute to gather — V' L with g^m 

L -(2+n-e )Qi =& a JZ n )(-V[)i n g Lm $ a AU + ^AU(-VL9L^aAU + *a,e(0^ n <,e(0(-^) 

=<AUU-vi 9Lm )<s> a AU + <s>aAU(~vig Lm )<P a AU + ^aAUi-vlg^AKAUin 

+ KAUL^BgL^aAU - ^aAUg^L^B. (A.9) 

We will track inverse powers of L, but otherwise group bounded operators together. We now introduce the 
following decomposition to commute —V^g^m through £ n localization. 

-VLql- =(-r^Vi)L- m L m ng Lm 
=(-r-X)£-™X T (x ro ) 2 

By Taylor's theorem, —r^ 1 V^ is an analytic function since Vl is, and V' L vanishes at r* = 0. For a localization 
$ being one of & a ,e, &b,e, or $ C:C , we write 

mU.-vig^AHHU.L^x^nfi-r-'vi)] 

=[a>(ej,x T (x„ l ) 2 ]i- m (-r; 1 y/) + r^ im [$(ej,-r; 1 y/. 9Lm ] 

=L m+n - 1 BL- m (-r- 1 V[) + r*g Lm L n - l B 
=L n ~ l B + r*L n - 1 B. 

By writing the decomposition as —gL m V[, the commutator can also be expressed as L n ~ 1 B + L n ^ 1 Br^,. 
Using this to continue with (|A.9[) gives 

L -(2+n-e) Ql =(-Vlg Lm )$ a AU(2KAZn)Zn + $«,e(0) + L^B + nL^B 

=(-VLg L ^cAU 2 + Ln ~' B + r*L n - x B 
=^cAU(-V[g L ^cAU + L n -'B + uL^B. 
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Taking the expectation value with respect to XaAi we get 

(x«V>,QiXaV>) =\\L^ (-Vlg Lm )H c , e (Ux a ipf + (x a i>,L 1+2n -*B Xa ip) + ( r * X a1>, L 1+2n ^B Xa ^) 
= ||L 2± ^(-^ 3im )^^(^) X ^|| 2 + 0(||L i ^ Xa VI| 2 ,i 1 ). 

The localizations can be commuted, and then the $> c ,e(£A localization can be replaced by e5<J> 0i£ (£ n ) local- 
ization. The same arguments as before show that all the commutators are lower order. 

>||£^*« ie (€J(-y£<7£»)Wll a + o(\\l^ x^W 2 ,^). 

Finally, we can commute the localizations back to get them in the standard orderer of our presentation. 
(XaV.QiXaV) ^L-^i-Vig^t^aM^X^ + O^L^x^f.L 1 ). 

We now consider the other terms in the expansion of the commutator with H3. All of these are lower 
order. To estimate the second term, Q2, we aim to absorb a factor of L rn in X|(x TO ). 

2i-( 2 +"-)Q 2 =-< !e (^[y i) ^ ra $ 0)e (0+^,e(0^<,e(0[^^„] 

This appears to be m powers of L too large to be lower order. We rewrite V' L = r*(r+ V£). 

2L -d+2n-e) Q2 = _ ^(^y^ir^vDxA^aAU + ^AU^xA^KAZMr^vL) 
= - KAU^x^AQaAUir^vl) + ZaAUxmXdxmXAZnXr^vL) 

- KAU^mX^AL^B + <f> a A€jL m X l (x m )L n - 1 B(r- 1 VL) (A.10) 
= - $ , 0>e ($a ) e(U)^„x m X i (x m )(r,- 1 ^) +x m X i (x m )$ , e (€ n )$ , 0fe (C n )(r- X) 

- B$ a . e (|„)L m+ "- 1 L"- 1 + $' 0je (^„)SL OT+ "- 1 + BL m+n ~ 1 . 

Since xX^(x) is bounded and m + n — 1 < 0, 

(x^Qax^) = (xaV>,£ 1+2 "^x^) = odlL^Wll 2 ^ 1 )- 

The third term is expanded using [[Vl,£„], £„] = V^(—L 2n ^ 2 ). We aim to absorb the derivative 9 r „ into 
the localization function e . 

L -(2+n-e)Q 3 = _ &>A£ n )Vl!L 2n - 2 d r ,g Lm $ a AU ~ ^MO^n^O^"" 2 

l-o»-«)q 3 = _ &>A$ n )d r ,vHg Lm <z> a AO - ^AUdL^KAUdr^ + B 

L -(l + 2n-e) Q3 = _ *£ e (£J£ n ^W*«,e(e„) - <&<M (€ n )fli» < )£ (€ n )€ n Yl + # 

Since £ (£) is bounded, Q3 = L 1+2n ~ e B, and Q3 is lower order in expectation value relative to Xai>- 
It is clear that Q4 is lower order without any rearrangement, 

2Q 4 =i 2+ "- t (-<, e (|J^'(-i 2 "- 2 ) 5 ^$ Q ,e(IJ + ^aAU9UKAUVl!(-L 2n - 2 )) 
=i 3n+m - e (< e (^)^ , X i ( Xm )$ a , £ (^)+$ a , £ (^)X i (x TO )< i£ (^)^)- 
The terms in Q5 are conjugates, up to a sign, so it is sufficient to estimate one of them. 

(</>, x^aAUL 2+n ^dr,gL^RzX^) =(V, x Q $a, e (^)i 2+ "- £ 9 rt5Lm i 3 ^ 3 i?XaV'>- (A.ll) 
Applying Cauchy-Schwartz and dropping bounded terms, we find this is bounded by 

ll^.Xa^H 2 + bL-i 4 "- Wll 2 £ ll^Xa^ll 2 + ||L 4n " W|| 2 = 0{\\L-^x a M\ 2 iL l ). (A.12) 
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The terms in Qg are dealt with similarly 

^,x a ®aAUL 2+n - e g' Lm R 3 x«^) =(V-,x Q $a, e (^)i 2+n - e i m ^i(x„ i )i 3n - 3 Sx^> 

= (XaV',i 4n+m - 1 - e i?XaV') 

^OiWL-^x^f.L 1 ). 
Gathering the result of the previous arguments, we have 

(V, [H,T n , m ]i>) >^aAUx^,L n ^{~d r J Lm d r , - ±g% m - C g L , n L 2 V[)<P a , £ (UXa fl 

This is very much like the term appearing in the Morawetz estimate, except for the scaling by L m in the 
function g^m and the factor of C on the potential term, and we proceed with a slight modification of that 
argument. Taking 

and x to denote characteristic functions, unrelated to Xa, of width 0(L~ m ) about r* = 0, and multiplying 
by L 3m , from an analogue of the argument leading to (|6.5p . we have 
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L 3m X v 2 d 3 n = - J L 3m r* X 'v 2 d 3 n + 2 J L 3m r* X vd r „vd 3 fi (A.13) 

< J L 5m rl X v 2 d 3 ii + J L m X (d rt v) 2 d 3 ^. (A.14) 
Continuing with the same argument, and multiplying the analogue of (|6.6p by L 2m , we have, 

L 2m v(0^) 2 d 2 Lu = - J L 2m X 'v 2 d 3 ii -2 J L 2m X vd r ,vd 3 V (A. 15) 

< J L 2m \ X '\v 2 d 3 f i + J L 3m X v 2 d 3 n + J L m x(d r ,v) 2 d 3 n)?( 

(A.16) 

Applying (|A.13[) to control the weighted L 2 norms, 

L 2m v(0) 2 <C ( J L 5m rl X v 2 dr, + J L m X (d rt v) 2 dr t ). (A.17) 
From integrating d rr g'[ m v 2 , the analogue of (|6.9p is 

As in ()6.10l) , the commutator can now be estimated 

aL 5m b 2 Corlxv 2 d 3 fi - G x f g Lm V[L 2 v 2 d 3 M . 



We now choose b sufficiently small, independently of L, so that for r* sufficiently small 

2 1 

- AabC Q >- 



(T + 1 £7 + 1 

-C^nVl >C <rb 2 r 2 
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We now use the assumption that m < 1/2. Since gL™ vanishes like r t L m , in the support of x, ~^9l™L 2 V[ > 
aL 5rn b 2 rlx, so that 

J 2g' Lm (d r ,v) 2 - l -g'£ m v 2 - Cg Lm V[L 2 v 2 d^ 

Since X^(x) and X±(x) are equivalent up to constant multiples, this gives the desired result. □ 

The previous lemma has localization in x m and £ n . To prove lemmas IA.1I and IA.21 we combine the 
localizations X|(x TO ) and X^{pc m ), in small and large x m , to get estimates localized in £ n . 

Proof of lemma \A.1[ From the estimate in the previous lemma, (|A.7j) , we have 

We divide gL m V[ into two parts based on whether x m is smaller or larger than one. Using \ to denote a 
smooth characteristic function on a scale of L~ m again, since —gL™L 2 Vl > L 5m r 2 x, by (|A.13|) . we control 
L 3m+n-€ - m tne re gi on Xm < i j n the region x m > 1, -g Lm L 2 V[ > —L 2 sgn(r*)V[ > -L 2 ~ m r* x V' L . Wc 
can combine both these estimates and take m = 1/2, 

(V, [H,r n>m }i>) >($ ,e(e„)XaV',i 3/2+ri - £ (-^- 1 ^)$ , e (e„)XaV'> 

The factor ~r~ x V' L is not really relevant, since it can be absorbed into Xa- To do this, we introduce a 
smooth, bounded function, /, which is equal to the inverse of {—r^Vf^) 1 ^ 2 on the support of Xon and then 
commute in the standard way. 

(iP,[H,T nim ]ip) >c\\L^ / ^(-r; 1 Vl)H a ^Ux a n+0(\\L^x a M\^L^ 
>c||L 1Z ^/(-r- 1 y/)3$ a , £ (| n )x^ll+0(||L i ^x a V'l| 2 ,i 1 ) 
>c||L 1Z ^$ aie (|J/(_ r -iy/)l Xc[ V'll+0(||L i ^ i x^l| 2 ,i 1 ) 
>c\\L ±I2 ^^ a Mn)x a n+0{\\L 1 - ± ^x a n 2 ,L 1 ) (A.18) 

Since §> a>e (x) > c$\ x \<i, we can replace the £ n localization to get 

[H,T n>m ]<P) >c\\L m ^^ in \< 1 x a n+0(\\L 1± ^x a n\L 1 ) 

□ 

To prove lemma IA.21 wc perform the same tricks of eliminating the x m localization to leave only £ n 
localization. This will require the introduction of a smooth approximation to ^ l- 5 <\x\<i- 

of lemma YA. e A To prove this result, we need a smooth approximation to ^ l- s <\x\<i- Let V^i be a smooth, 
weakly increasing function with support on [1/2, oo) and identically one on [l,oo), and let '5 2 be a smooth, 
weakly decreasing function with support on [0,2] and identically one on [0,1]. We extend these as even 
functions. Let 

*(A- 5 ,0 = *(A-^)*2(£)- 
We first show ^ dominates ^ \-s<\»\<i- Decomposing to harmonics and applying the spectral theorem 

in 

>X[a-s,oo)(z)X[o,i](z) 

>*A-«<|x|<i, (A.19) 
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so that for ip g L 2 (97l), 
We also note 



ii*(^,a^ii 2 >ii*L-<i«j<i^n 2 . 



\- 5 *(\- S , < £$ c , £ (£) < £4>a,e(£). (A.20) 

Since X P(L _<5 ,^„) = V&i(^ rl+ ^) v I / 2(^„), by the order reduction lemma, following the same argument which 
justifies equation (|A.4[) . if F x is a function with bounded derivative, such as Xy{x), X|-(x), or (— r~ 1 V[), or 
a smooth compactly supported function, 

[F 1 (x TO ),*(i- 5 ,^)]=[F 1 (x TO ),* 1 (^ +5 )]* 2 (C n )+*i(^ +a )[F 1 (x m ),* 2 (^)] 
=L m+,i+ ^ 1 S* 2 (|„) + * 1 (^, i+5 )L m+ "" 1 B 

We start with the estimate involving the derivative operators. We relax the estimate by a factor of 
JTj(x„)2 and write d r , in terms of £ n before beginning to commute the localization, which generates the 
standard lower order terms. 

(V>, [H 1+ H 3 ,T n ^) >||L™-tl i (x„)^ rt $ a£ (^) XQ ^l| 2 + 0(||L i ^ X ^|| 2 ,L 1 ) 

^H^-^KXnKn^^JXoV'll 2 + OiWL^X^W 2 , L 1 ) 
>||L 1 -*|„$ Q , e (|„)X 1 (x„)XaV'l| 2 + 0(||i^Xa^|| 2 , L l ). 

We now use equation (|A.20[) to introduce ^(X^ 5 , £„) and continue commuting localizations. 

+ £r 3 ,r n ,„]v) ^Ili^^-^^-^uxxCx^x^^lp + odlL^x^ll 2 ^ 1 ) 

^-^X^^iL-^UXc.ir + OillL^x^^L 1 ) (A.21) 

The last commutator is of order 2n + 8 — 1 in L. Since the introduction of 4" reduced the order of the entire 
expression to 1 — e/2 — 5, the term from the commutator is still lower order, of order 2n < 1/2 + n. 
Turning to the potential term in equation (|A.7|) . we have 

(V>,[ffi+ff 3 ,r„,„]v.>} >||L 1 -4x T (x n )(-r7 1 y/)^ a , £ (|„)xa^|| 2 + o(||i i ^ i XaV'll 2 ,i 1 )- 

Using the same argument as in (|A.18|) . we can eliminate the factor of (— r* 1 V[)^. 

(V, [Hi + H 3 ,T n ^j) >||L 1 -tX T (x n )<i> a:e (OxaV'l| 2 + OiWL^x^Al 2 , L 1 ). 

To continue, we drop a factor of L s and then start to commute localizations, and use <J> ai£ (x) > 4'(A~' 5 , x). 
We need to drop the factor of 5 powers of L to control commutators involving ^(L~ s ,£„). 

[H 1+ H^Y n , n ]iP) >||L 1 -*- a JT T (x n )* o , e (€ Tl )x«^l| a + 0(||i i¥a Xa^|| 2 ,i 1 ) 

>||L 1 -t- 5 $ a ,e(^JX T (x n )XaV'l| 2 + C»(l|i i ^ 2 XaV'l| 2 I i 1 ) 
>||i 1 -f-^(X- 5 ^ n )X T (x n )XaVl| 2 + 0(||i 1 ^XaV'l| 2 1 i 1 ) 

^ll^-i-'XTCxn^CL-^^Jx^f+OaiL^Xa^ll'.i 1 ). (A.22) 
We use ()A.2p to combine (|A.21|) and (|A.22j) . and use (|A.6|) to control the contribution from iJ 2 . 

($,[H,T n , n }<V) >||i 1 -*-^(i- 5 ,OxaV'll 2 + 0(||i i ^ i XaV'll 2 ,i 1 )- 

We conclude by using (|A.19[) to replace ^(L~ s ,^ n ) by |<i, 

^,[H,T n<n ty) >||L 1 -^^ i -,<| €J < 1 x Q V!| 2 + 0(||L^ 1 x Q V'l| 2 ,i 1 )- 

□ 



38 



Acknowledgement 

We would like to thank J. Sterbenz for providing valuable discussions concerning the results in [6]. 

References 

[1] A. Bachelot, J. -P. Nicolas, Equation non lineaire de Klein-Gordon dans des metriques de type 
Schwarzschild. C. R. Acad. Sci. Paris Sr. I Math. 316 (1993), no. 10, 1047-1050. 

[2] P.H. Berard, Spectral Geometry: Direct and inverse problems. Lecture notes in Mathematics, Springer 
Verlag 1207 1986. 

[3] P. Blue and A. Softer, Semilinear wave equations on the Schwarzschild manifold I: local decay estimates, 
Adv in Differential Equations, 8 (2003) 595-614. 

[4] P. Blue and A. Softer, The wave equation on the Schwarzschild metric II: local decay for the spin 2 
Regge Wheeler equation, J. Math. Phys. 46 (2005), 012502, 9 pp. 

[5] P. Blue and A. Softer, Phase Space Analysis on some Black Hole Manifolds, math. AP/051 1281, preprint. 

[6] P. Blue and J. Sterbenz, Uniform Decay of Local Energy and the semilinear wave equation on 
Schwarzchild space, Commun. Math. Phys., to appear. 

[7] N. Burq, P. Gerard, N. Tzvctkov, Strichartz inequalities and the nonlinear Schrddinger equation on 
compact manifolds. Amer. J. Math. 126 (2004), no. 3, 569-605. 

[8] T. Cazenave, Semilinear Schrddinger equations. Courant Lecture Notes in Mathematics, 10. New York 
University, Courant Institute of Mathematical Sciences, New York; American Mathematical Society, 
Providence, PJ, 2003. 

[9] D. Christodoulou and S. Klainerman, The global nonlinear stability of the Minkowski space. Princeton 
Mathematical Series, 41. Princeton University Press, Princeton, NJ, 1993. 

[10] M. Dafermos, The interior of charged black holes and the problem of uniqueness in general relativity. 
Comm. Pure Appl. Math. 58 (2005), no. 4, 445-504. 

[11] M. Dafermos and I. Rodnianski, The red-shift effect and radiation decay on black hole space times, 
gr-qc/0512119, preprint. 

[12] J. Dimock and B.S. Kay, Classical and quantum scattering theory for linear scalar fields on the 
Schwarzschild metric. I. Ann. Physics 175 (1987), no. 2, 366-426. 

[13] S.-I. Doi, Smoothing effects for Schrddinger evolution equation and global behavior of geodesic flow. 
Math. Ann. 318 (2000), no. 2, 355-389. 

[14] G.F.R. Ellis and S.W. Hawking Ellis, The large scale structure of space-time. Cambridge Monographs 
on Mathematical Physics, No. 1. Cambridge University Press, London-New York, 1973. 

[15] G.F.R. Ellis and K.S. Virbhadra, Schwarzschild black hole lensing, Phys. Rev. D. 62 (2000) 084003 

[16] L.C. Evans, Partial differential equations. Graduate Studies in Mathematics, 19. American Mathematical 
Society, Providence, RI, 1998. 

[17] F. Finster, N. Kamran, J. Smoller, S.-T. Yau, Decay of Solutions of the Wave Equation in the Kerr 
Geometry, Commun. Math. Phys. 264 (2006) 465-503. 

[18] J. Ginibre and G. Velo, Conformal invariance and time decay for non linear wave equations. I. Ann. 
Inst. Henri Poincare, 47 (1973) 221-261. 

[19] D. Hairier, Completude asymptotique pour I'equation des ondes dans une classe d'espaces-temps station- 
naires et asymptotiquement plats. Ann. Inst. Fourier (Grenoble) 51 (2001), no. 3, 779-833. 



39 



[20] D. Harrier and J.-P. Nicolas, Scattering of massless Dirac fields by a Kerr black hole. Rev. Math. Phys. 
16 (2004), no. 1, 29-123. 

[21] A. Hassell, T. Tao, and J. Wunsch, A Strichartz inequality for the Schrdinger equation on nontrapping 
asymptotically conic manifolds. Comm. Partial Differential Equations 30 (2005), no. 1-3, 157-205. 

[22] I. Laba and A. Soffcr, Global existence and scattering for the nonlinear Schrddinger equation on 
Schwarzschild manifolds. Helv. Phys. Acta 72 (1999), no. 4, 274-294. 

[23] R. B. Lavine, Commutators and scattering theory. I. Repulsive interactions. Comm. Math. Phys. 20 
(1971), 301-323. 

[24] H. Lindblad and I. Rodnianski, Global existence for the Einstein vacuum equations in wave coordinates. 
Comm. Math. Phys. 256 (2005), no. 1, 43-110. 

[25] C. Morawetz, Time decay for the nonlinear Klein-Gordon equations. Proc. Roy. Soc. Ser. A 306 (1968) 
291-296. 

[26] J.V. Ralston, Solutions of the wave equation with localized energy, Comm. Pure. Appl. Math., 22 (1969), 
807-823. 

[27] T. Regge and J. A. Wheeler, Stability of a Schwarzschild singularity, Phys. Rev. 108 (1957), 1063-1069. 

[28] H.F. Smith and CD. Sogge, Global Strichartz estimates for nontrapping perturbations of the Laplacian. 
Comm. Partial Differential Equations 25 (2000), no. 11-12, 2171-2183. 

[29] A. Soffer, On the many-body problem in quantum mechanics. Astrisquc 207, (1992), 6, 109-152. 

[30] CD. Sogge, Lectures on nonlinear wave equations, Monographs in Analysis, II. International Press, 
Boston, MA, 1995. 

[31] W.A. Strauss, Nonlinear wave equations. CBMS Regional Conference Series in Mathematics, 73. AMS, 
Providence, RI, 1989. 

[32] M.E. Taylor, Tools for PDE. Pseudodifferential operators, paradifferential operators, and layer poten- 
tials. Mathematical Surveys and Monographs, 81. American Mathematical Society, Providence, RI, 
2000. 

[33] S.A. Tcukolsky, Rotating black holes: Separable wave equations for gravitational and electromagnetic 
perturbations. Phys. Rev. Lett. 29 (1972), 1114-1118. 

[34] B.F. Whiting, Mode stability of the Kerr black hole. J. Math. Phys. 30 (1989), no. 6, 1301-1305. 

[35] F.J. Zcrilli, Effective potential for even-parity Regge- Wheeler gravitational perturbation equations. Phy. 
Rev. Lett. 24 (1970) 13, 737-738. 



40 



